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ABSTRACT 


I 

An  expemded  package  of  interactive  FORTRAN  con^uter  pro- 
greuos  has  been  developed  for  the  emalysis  and  forecasting 
of  time  series  data.  The  package,  called  the  Time  Series 
Editor,  is  designed  to  employ  the  iterative  Box- Jenkins 
methodology  of  time  series  emalysis.  The  Time  Series  Editor 
was  developed  for  time-shared  use  on  the  Control  Program/ 

Cambridge  Monitor  System  (CP/CMS)  at  the  U.S.  Naval  Post- 
gradviate  School,  but  can  be  modified  for  use  on  other  time- 
sharing systems  with  a FORTRAN  capediility.  The  Time  Series 
Editor  assists  in  data  preparation  and  entry,  euialysis, 
modeling,  forecasting  and  diagnostic  testing.  Utilization 
of  the  package,  following  the  included  User's  Guide,  requires 
only  a limited  knowledge  of  the  computer  system,  with  all 
required  user  responses  interactively  prompted  by  the  Editor. 
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Operations  researchers,  military  plemners  and  programmers, 
statisticiems , economists,  marketing  personnel,  managers  emd 
others  are  often  faced  with  the  need  to  analyze  data  in  the 
form  of  a time  series,  which  can  be  thought  of  as  a sequence 
of  observations  of  either  a deterministic  or  stochastic 
process.  In  most  cases,  the  objective  of  the  analysis  is 
to  determine  patterns  or  other  recognizable  behavior  appar- 
ent in  the  data,  and  to  then  formulate  a suitable  mathemati- 
cal model  for  the  time  series  from  which  forecasts  of  futxire 
behavior  can  be  obtained.  The  value  to  a decision  maker 
of  being  able  to  predict  the  future  with  some  reasonable 
and  statistically  qucmtifi2d}le  degree  of  accuracy  cem  not 
be  overstated.  For  example,  in  areas  such  as  budget  expen- 
ditures, recruiting  performance,  commodity  prices,  population 
levels,  resource  consumption,  msuipower  levels  2md  consumer 
demand  for  products,  decision  maUcers  chaurged  with  plauining 
for  the  future  should  base  their  decisions  at  least  in 
part  on  the  best  availad^le  predictions  about  the  futvire 
behavior  of  the  time  series  in  question. 

Until  the  late  1960 's,  the  techniques  enployed  in  time 
series  analysis  were  primarily  those  of  spectral  analysis, 
with  applications  of  haunnonic  analysis  and  mathematical 
transform  theory.  For  an  introduction  to  these  topics,  the 
reader  is  referred  to  Anderson  [Ref.  2] , Lewis  [Ref.  7] 
amd  Jenkins  [Ref.  10].  Due  to  the  high  degree  of  mathematical 
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sophistication  required  by  the  spectral  analysis  approach, 
the  capedjility  to  perform  time  series  analysis  resided  nearly 
exclusively  with  mathematiflians  and  electrical  engineers. 

As  a consequence,  the  majority  of  decision  meUcers  came  to 
use  more  understandable  (euid  far  less  powerful)  methods 
such  as  simple  moving  averages  or  exponential  smoothing. 
However,  in  the  late  1960's  and  early  1970* s,  the  statistical 
analysis  of  time  series  by  the  methods  developed  by  Box 
and  Jenkins  [Ref.  4]  has  gained  widespread  acceptance. 
Although  these  methods  are  mathematically  nearly  equivalent 
to  the  spectral  approach  (there  are  transformations  that 
interconnect  the  two  methodologies) , the  Box-Jenkins  approach 
is  described  in  vocabulary  more  fcuniliar  to  operations 
researchers,  statisticians,  economists  and  managers,  and 
is  therefore  being  used  more  and  more  by  them  in  building 
models  from  measurements  of  the  past. 

Many  algorithms  and  computer  programs  for  performing 
the  analyses  required  by  the  Box-Jenkins  approach  have  been 
developed  and  are  available  from  several  sources.  Perhaps 
the  best  source  is  the  collection  of  FORTRAN  computer  sub- 
routines which  resides  in  the  International  Mathematical 
and  Statistical  Library  (IMSL)  [Ref.  9].  The  primary  prob- 
lem with  using  the  available  computer  resources  lies  not 
in  any  deficiency  of  the  programs  themselves,  but  with  the 
veiry  nature  of  the  Box-Jenkins  methodology  itself.  The 
Box-Jenkins  method  is  an  iterative  approach  to  time  series 
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analysis,  which  is  briefly  described  in  Figure  1.  [See 
Wheelwright  and  MakrideUcis,  Ref.  16.] 


As  indicated  by  Figure  1,  the  Box- Jenkins  method  is  a 
multi-stage,  iterative  process.  It  begins  with  the  postu- 
lation of  a general  class  of  time  series  models  which  has 


been  found,  experimentally,  to  be  extremely  rich.  There- 
after, the  procedure  continues  as  a trial-and-error  process, 
with  decisions  points  where  the  analyst  is  reqiiired  to 
select  the  next  direction  based  on  the  best  information 
available  to  him.  Since  each  stage  of  the  process  outlined 
in  Figure  1 may  consist  of  several  sub-steps,  tjie  modeling 
process  itself  can  become  quite  time-consuming,  even  with 
the  ready  availability  of  the  IMSL  softwaure  resources.  For 
example,  an  analyst  employing  the  IMSL  subroutines  in  a 
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FIGURE  1.  The  Box-Jenkins  Forecasting  Methodology 
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5.  Write  a progreun  to  call  the  IMSL  subroutine  which 
transforms  the  time  series  to  adjust  for  seasonal 
patterns,  nonstationaury  behavior  or  other  behavior 
which  deviates  from  that  assumed  for  the  class  of 
models  postulated. 

6.  Repeat  steps  (2)  through  (4)  with  the  tramsformed 
data. 

7.  Review  the  statistical  properties  of  the  auto- 
correlations euid  partial  autocorrelations  for 
tentative  model  identification. 

8.  Write  a program  to  call  the  IMSL  program  that 
estimates  the  model  parameters  and  computes  the 
residuals. 

9.  Write  a program  to  perform  goodness-of-fit  tests 
for  the  model. 

10.  Analyze  the  model  residuals,  using  steps  (1)  through 
(9) , as  with  the  original  time  series. 

11.  Refine  the  model  using  information  gained  by 
examination  for  any  structure  in  the  residuals. 

12.  Repeat  the  preceding  process  until  no  structure 
remains  in  the  residuals. 

13.  When  an  adequate  model  has  been  obtained,  write  a 
program  to  call  the  IMSL  subroutine  which  forecasts 
and  determines  confidence  intervals  for  future 
values  of  the  time  series. 

Between  each  successive  pair  of  steps,  the  emalyst  must 
numually  intervene  in  the  analysis  process,  meUcing  a decision 
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subjectively  based  on  the  information  then  available. 

Thus,  much  analyst  interaction  is  necessary  in  order  to 
determine  a suitable  mathematical  model  and  forecast  equa- 
tion. Even  with  rapid  computer  job  turnarotmd  time,  the 
process  described  above  could  easily  consume  a full  working 
day  or  more. 

This  paper  describes  an  effort  to  alleviate  some  of  the 
organizational  problems  inherent  in  the  modeling  of  time 
series  using  the  Box- Jenkins  iterative  approach.  This 
effort  consists  of  an  interactive  con^uter  program  package 
which  provides  access  in  a structured  way  to  the  most  useful 
IMSL  and  other  subroutines  necessary  for  Box-Jenkins  time 
series  analysis.  This  package,  called  the  Time  Series  Editor, 
is  written  for  use  with  the  Naval  Postgraduate  School's 
Control  Progam/Cambridge  Monitor  System  (CP/ CMS) . Since  all 
working  programs  with  the  exception  of  the  executive  routine 
itself  are  written  in  FORTRAN,  the  Time  Series  Editor  should 
be  adaptable  for  use  on  other  PORTRAN-capable  time-sharing 
systems.  The  Time  Series  Editor  assists  the  analyst  in 
data  prepauration  and  entry,  model  construction  amd  diag- 
nostic testing,  and  time  series  forecasting.  In  fact, 
with  the  User's  Guide  provided  as  an  Appendix  to  this  report, 
a complete  Box-Jenkins  time  series  analysis  can  be  per- 
formed in  a short  time  by  even  a no'«ice  computer  user. 

For  this  reason,  the  Time  Series  Editor  could  be  a valuable 
instructional  tool  for  laboratory  use  in  a time  series 
course. 
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A brief  overivew  of  Box-Jenkins  methodology  is  provided 
in  Chapter  II,  to  serve  as  a point  of  reference  for  the 
material  which  follows.  Chapter  III  contains  a description 
of  the  algorithm  employed  to  estimate  non-lineau:  least- 
squares  parameters  for  a generalized-  (seasonal,  nonseasonal, 
stationary,  or  nonstationary)  Box-Jenkins  model.  Chapter  IV  | 

contains  descriptions  of  each  of  the  major  programs  in  < 

the  Time  Series  Editor.  Chapters  V and  VI  conted.n  examples  | 

of  actual  use  of  the  Time  Series  Editor,  using  a non-  | 

seasonal  and  a seasonal  series,  respectively.  Chapter  VII  ^ 

includes  a summary  of  the  report,  and  recommendations  for  | 

future  additions  to  the  Time  Series  Editor.  The  Appendices  | 

include  a User's  Guide  to  the  Time  Series  Editor,  with  j 

1 

instructions,  san^Jle  user  sessions  and  sample  outputs,  as  | 

well  as  listings  of  the  Editor  itself  and  all  primeary  | 

subprogrcun  packages . | 
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I I . BOX- JENKINS  METHODOLOGY:  A PRIMER 


This  chapter  provides  a brief  overview  of  the  methodology 


developed  for  the  2malysis  of  time  series  data  by  G.E.P. 


Box  eind  G.M.  Jenkins,  in  order  to  facilitate  understanding 


of  the  models  and  programs  en^loyed  in  the  Time  Series 


Editor.  For  more  detail  on  the  material  covered  in  this 


chapter,  the  interested  reader  is  referred  to  the  texts  by 


Anderson  [Ref.  1],  Box  amd  Jenkins  [Ref.  4],  Pindyck  and 


Rubinfeld  [Ref.  15],  Nelson  [Ref.  14],  and  Mabert  [Ref.  12]. 


A.  LISTING  OF  NOTATION 


The  following  is  a listing  of  the  notation  employed 
in  this  chapter  and  in  the  remainder  of  this  report.  Each 


symbol  is  provided  with  a brief  title  or  explanation  of  its 


a time  series 


the  length  of  the  time  series  (number  of  terms) 


estimates  of  the  meam  auid  variamce  of  the  time 


series  {z^} 


transform  of  {z^},  where  z^' 


differenced  series  from  {z  },  where 


autocorrelation  of  lag  k 


estimate  of  autocorrelation  at  lag  k 


backward  difference  operator  of  order  d 


•9.  ■ 


s 

«(> 


e 

0 


S($,0) 


p 


p 


q 


Q 

'Kb) 


seasonal  backward  difference  operator  of  order 
D with  season  s 

backward  shift  operator  of  order  s 

length  of  the  time  series'  season 

non- seasonal ^ autor egr es s i ve  parameters 
(estimates,  $) 

seasonal  autoregressive  parameters  (estimates,  $) 

non-seasonal^moving  average  parameters 
(estimates,  8) 

A 

seasonal  moving  average  parameters  (estimates,  9) 
series  of  white  noise  terms  of  the  process 
variance  of  the  white  noise  process 

I 

sum  of  squared  residuals  of  the  model 

number  of  non-seasonal  autoregressive  parameters 

number  of  seasonal  autoregressive  parameters 

nvimber  of  non-seasonal  moving  average  parameters 

number  of  seasonal  moving  average  parameters 

transfer  function  of  a linear  filter 

partial  autocorrelation  of  log  k 


B.  TIME  SERIES,  STOCHASTIC  MODELS  AND  STATIONARITY 

A distinguishing  feature  of  modem  time  series  amalysis 
is  that  the  sequence  of  observations  of  a given  varieUsle  is 
considered  to  be  a realization  of  jointly  distributed 


I 


I 


I 
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which  defines  the  level  about  which  it  fluctuates  in  tiioe, 
and  a constant  variance 


Ef  (2^  - Ji)  1 


/ (z  - ) p(z)dz  , 


which  measures  its  variability  eibout  the  level  of  the  mean. 
Since  the  probability  distribution  p(z)  is  the  seuae  for  all 
times  t,  its  form  cem  be  discovered  by  plotting  a histogram 
or  relative  frequency  plot  of  the  set  of  observations 
(z^,Z2r • . . . Additionally/  the  mean  and  variance  of 
the  stochastic  process  can  be  estimated  by  the  sample 
averages  taken  over  time: 


I I ^ 


5 I - 


C.  THE  BASICS  OF  DIFFERENCE  OPERATORS  AND  LINEAR  FILTER 
MODELS 

This  section  provides  a brief  discussion  of  the  differ- 
ence operators  currently  employed  in  writing  Box-Jenkins 
models.  Those  most  commonly  used  are  described  below  in 
Table  I. 


I 


Operator 


Name 


Definition 


B 


m 


Baclcweurd  Shift  Operator  Bz  * z . 
of  order  m ^ 


■ *t-s 


Forweurd  Shift  Operator 
of  order  m 


F - B 


-1 


''■'t  * 'tw 


S 


Backward  Difference 
Operator  of  order  D 
and  season  s 


7°  - (1  - B)° 


7°  - (1  - B®)° 


7z^  - zt  - 2^-1 


Ted}le  I.  Basic  difference  operators. 


It  is  noted  here  that  7 z^  could  be  written  as  a 


binomial  exp2uision;  that  is. 


7°z. 


(1  - B)‘'z, 


( I (j)  {1)^(-B)°“^)z^ 
j-0 


I (^)(-l)°‘^Zt-D+j 
j-0 


The  stochastic  models  employed  in  the  Box-Jenkins  method 
are  based  on  the  idea  that  a time  series  in  which  successive 


2 


I 

values  are  highly  dependent  can  be  regarded  as  having  been 

generated  from  a series  of  independent  "shocks”  a^.  These 

shock  terms,  sometimes  referred  to  as  white  noise  terms, 

are  random  drawings  from  some  fixed  distribution,  generally 

2 

assumed  Normal  with  meeui  zero  emd  variance  a . The  white 

ci 

noise  process  a^  is  supposed  trsmsformed  to  the  process  z^ 
by  what  is  known  as  a linear  filter,  shown  in  Figure  2 below 
[Ref.  4,  p.  8] . The  linear  filtering  operation  takes  a 
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■ 1 + ••• 

is  the  linear  operator  that  transforms  a^  into  z^,  called 
the  transfer  function  of  the  filter.  The  sequence 

formed  by  the  weights  may  be  (theoretically) 
finite  or  infinite;  if  this  sequence  is  finite,  or  infinite 
2uid  converges,  the  filter  is  said  to  be  steUble  and  the 
process  z^  stationary.  The  peirameter  p is  then  the  value 
about  which  the  process  varies.  Otherwise,  z^  is  called 
nonstationeury,  and  y has  no  specific  meeuiing  except  as  a 
reference  point  for  the  level  of  the  process. 

D.  HOMOGENEOUS  NON-STATIONARY  AND  SEASONAL  PROCESSES 

In  m2uiy  cases  of  interest  to  the  time  series  analyst, 
the  process  being  exeunined  is  not  stationary.  Instead, 
the  probaUsilistic  structure  of  the  process  generating  the 
time  series  changes  with  time.  For  example,  there  may  be 
a trend  or  seasonal  pattern  in  the  time  series.  However,  if 
the  series  can  be  found  to  exhibit  behavior  which  is  somewhat 
homogeneous,  then  the  series  can  frequently  be  transformed 
into  a stationeury  series  that  can  be  described  by  Box-Jenkins 
models . 

A series  is  said  to  be  homogeneous  nonstationary  of 
order  d if  w^  * ^ stationary  series  for  some  integer 

d.  Here  7 = (1  - B)  denotes  the  difference  operator  such 
that  7^a^  » z^  - The  fact  that  a given  series  is 

nonstation^u:y  can  be  determined  quickly  by  exzurdnation  of 
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the  autocorrelation  function,  to  be  discussed  later.  The 
series  to  be  modeled  should  be  differenced  until  the  re- 
sulting series  appears  stationeury  or  until  the  procedure 
appears  to  be  m2dcing  no  improvement. 

A time  series  is  called  seasonal  when  it  exhibits 
cyclical  behavior  over  time,  showing  a regular  period. 

Exan^les  of  seasonal  processes  include  monthly  rainfall, 
monthly  crop  yields,  sale  volume  of  bathing  suits,  livestock 
production  rates  and  energy  consumption.  Seasonal  patterns 
are  often  easy  to  spot  simply  by  observing  a plot  of  the 
series,  or  through  knowledge  of  the  process  that  generated 
the  series.  However,  in  many  cases  where  the  variadiility  of 
the  series  is  large,  seasonal  patterns  will  be  difficult  to 
distinguish  from  other  fluctuations.  Recognition  of ‘season- 
ality in  a series  is  important,  since  it  provides  information 
useful  in  modeling  auid  forecasting;  as  will  be  shown,  the 
autocorrelation  ftmction  can  assist  in  recognizing  seasonality 
The  Box-Jenkins  modeling  approach  for  seasonal  nonstationary 
time  series  of  season  length  s is  to  first  transform  the 
series  using  eui  appropriate  seasonal  differencing  operator 
to  "sweep  out"  the  seasonal  effect;  that  is,  take 


If  there  is  a trend  or  other  types  of  non-stationarity  pre 
sent,  as  well  as  the  seasonal  effect,  the  series  can  be 
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differenced  again  with  the  ordinary  difference  operator 
until  stationarity  is  achieved.  This  treuisform  would  be 
written 


w =<  V°  z 
t St 

The  transformed  series  is  then  modeled  as  a stationaury 
series.  For  seasonal  series  where  con^lete  stationarity 
cem  not  be  achieved,  a class  of  seasonal  models  is  avail-' 
aUsle;  these  are  discussed  later  in  this  Chapter. 

E.  GENERAL  CLASSES  OF  BOX- JENKINS  MODELS 

This  section  will  provide  brief  descriptions  of  several 
of  the  most  commonly  used  classes  of  Box-Jenkins  models, 
including  autoregressive  (AR)  models,  moving  average  (MA) 
models,  mixed  and  integrated  autoregressive-moving  average 
(ARMA/ARIMA)  models,  and  a generalized  seasonal  ARIMA  model. 
1.  Autoregressive  Models 

In  autoregressive  (AR)  models,  the  current  value  of 
the  series  under  study  is  expressed  as  a linear  combination 
of  previous  series  values  that  explain  the  current  obser- 
vation, plus  an  vinexplained  random  (white  noise)  term  a.  . 

For  a stationary  series,  or  one  that  has  been  treuisformed 
to  stationarity,  the  deviation  from  the  mean  value  for  each 
period,  =*  ^t  ” modeled  as  dependent  on  weighted 

values  of  the  previous  deviations  from  the  meem.  The  AR 
model  C2U1  be  written 


m 


+ ...  + <fr  z . 

P t-p 


In  the  expressions  which  follow  the  tilde  will  be  dropped, 
but  it  is  to  be  understood  that  the  mecui  has  already  been 
subtracted  from  each  observation.  This  equation  is  called 
an  autoregressive  process  of  order  p,  where  (j>j  represents 
a scalar  weighting  coefficient  for  the  previous  period. 
The  model  ceui  also  be  written  in  more  compact  form  as 


(J)(B)z^  * a^  , 


where 


({)(B)  = (1  - (}.J^B  - - ...  - 


a polynomial  in  the  operator  B.  Then,  (B)  represents  a 
general  non-seasonal  AR  operator.  The  reason  these  models 
are  called  autoregressive  can  be  seen  from  recalling 
regression  emalysis;  here  the  model  relates  the  dependent 
variable  z^  to  a set  of  explamatory  variedsles 
...,z^_p),  which  are  previous  values  of  the  time  series. 
Therefore,  the  model  is  autoregressive.  For  ex2unple,  when 
the  value  of  the  time  series  at  time  t depends  only  on  the 
value  at  time  t-1,  the  model  becomes  an  autoregressive  model 
of  order  1,  denoted  AR(1).  Similarly,  if  the  value  in 
period  t depends  on  the  values  in  periods  t-1,  t-2  2md  t-3. 


In  moving  average  (MA)  models,  the  assumption  is  made 


that  the  current  value  of  the  time  series  at  period  t can 
be  expressed  as  a linear  combination  of  the  previous  fore- 
cast errors  (or,  residuals),  a^  “ ^t  ” *t  * general 

equation  for  this  model  can  be  written 


2.  — a.  — ® 1 ® O ""  • • • — ^ 

t t 1 t-1  2 t-2  q t-q 

which  is  called  a moving  average  model  of  order  q,  vdiere 

ttl 

represents  the  scalar  weighting  coefficient  for  the  q 
previous  period.  Like  the  AR  model,  the  MA  model  caui  be 
written  in  more  con5)act  operator  notation  as  follows: 

2^  = 0(B)  a^  , 


a polynomial  in  the  operator  B.  The  MA  model  implies  that 
the  emalyst  can  gain  valuable  information  for  future 
predictions  by  consideration  of  the  weighted  sum  of  a 
number  of  previous  forecast  errors  or  residuals.  The 
following  eure  examples  of  MA  models  of  order  1 2uid  order 
3,  MA(1)  and  MA(3) : 


MA(1)  : - (1  - 0j^B)a^  = a^  - ? 


MAt3)  : Z^  » (l-6j^B-e2B^-03B^)a^ 


*t  “ ®l*t-l  ” ®2*^t-2  “ ®3^t-3 


3.  Mixed  and  Integrated  Autoregressive-Moving  Average 
Models 

A natural  extension  of  the  autoregressive  and  moving 
average  models  is  to  construct  a combination  of  the  two. 

Such  mixed  processes  cure  called  autoregressive-moving  average 
(ARMA)  processes,  of  orders  p and  q,  often  written  as 
ARMA(p,q).  Box  and  Jenkins  [Ref.  4]  notes  that  for  many 
series  encountered  in  practice,  the  inclusion  of  both  AR 
and  MA  terms  in  a model  results  in  fewer  total  parameters 
than  would  be  required  for  either  a pure  AR  or  pure  MA 
process.  The  ARMA(p,q)  model  may  be  written  in  the  usual 
form  as 


•^St-l  '*’2*t-2  “^p^t-p  S " ®lS-l  " ®2^t-2  “ • * • 


-L-LIS^JUUK 


in  the  more  contact  operator  notation/  the  same  model  is 
written 


4)(B)z^  * 0(B)  . 

As  an  exan^le,  eui  ASMA.(2,2)  model  would  be  written  as: 

(1- <|>j^B- ((i2B^)z^  » (1- 0j^B- 02B^)a^  / 

or  expemded/  as: 

z^  » <l»i2^  1 + - 0,a^  , - 0-,a^  - . 

t 1 2 't"*2  u X t“*X  2 t*2 

Recall  that  with  differencing  or  other  transformation 
techniques/  a stationary  series  can  be  obtained  from  a 
homogeneous  (seasonal  or  non-seasonal)  non-stationeary  series. 
When  the  differencing  technique  is  integrated  directly  into 
the  ARMA  model/  the  result  is  an  integrated  autoregressive- 
moving  average  (ARIMA)  model.  The  model  is  often  written 
ARINA(p/d/q) / where  p amd  q retain  their  previous  meamingS/ 
amd  d represents  the  order  of  (non-seasonal)  differencing. 

The  model  may  also  be  generalized  by  including  an  additional 
constant  tem  0^  that  will  have  the  effect  of  adding  a 
deterministic  polynomial  trend  of  order  d.  The  value  of  0^ 
may  be  found  in  general  by  the  formula 
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The  "differenced"  form,  equation  (1) , is  the  form  usually 
seen.  In  the  differenced  form,  the  transfer  function  <|>(B) 
is  assumed  to  be  stable;  that  is,  all  the  roots  of  the 
difference  equation  (^(B)  » 0 are  outside  the  unit  circle. 
Clearly,  ^(B)  is  not  a transfer  function  of  a stationary 
series  when  d > 0,  since  f (B)  has  d roots  on  the  unit 
circle. 
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The  general  class  of  ARIMA  models  previously  dis- 
cussed can  be  further  generalized  to  allow  for  the  modeling 
of  seasonal  models  with  the  addition  of  appropriate  seasonal 


nonseasonal  and  seasonal  operators 


F.  THE  AUTCX:ORRELATION  AND  PARTIAL  AUTOCORRELATION  FUNCTIONS 


The  autocorrelation  function  of  a time  series  describes  the 
association  (mutual  dependence)  cunong  values  of  the  seune 


1 


series  taJcen  at  different  time  periods,  the  difference  in 
time  being  referred  to  as  the  lag,  denoted  k.  Autocorrela- 
tion of  itself  in^lies  nothing  edxjut  a change  in  one 
veuriable  causing  a change  in  2uiother.  However,  the  auto- 
correlations may  provide  in^ortant  intormation  2dx>ut  the 
structiire  of  a data  set  and  its  pattern.  In  a set  of  com- 
pletely rauidom  data  the  theoretical  autocorrelations  among 
successive  values  will  always  be  zero,  whereas  data  values 
of  strong  seasonal  or  cyclical  behavior  will  be  highly 
autocorrelated.  This  can  be  estimated  using  the  time- 
averaged  sample  autocorrelation,  denoted  r^^,  defined  as: 


N-k 


N t»l  ^ 


- z)  (z 


t+k 


- z) 


A plot  of  the  autocorrelation  function  versus  the  lag  k, 
called  a correlogreun,  is  very  useful  for  the  purpose  of 
determining  if  a process  is  stationary  and  for  identifying 
the  appropriate  model. 

Another  function  useful  in  identifying  the  appropriate 
model  for  a given  series  is  the  partial  autocorrelation 
function,  often  called  the  pauto  function.  Partial  auto- 
correlations are  analogous  to  autocorrelations  in  that  they 
indicate  the  relationship  of  the  values  of  a time  series  to 
various  time-lagged  values  of  the  same  series.  However, 
they  differ  from  autocorrelations  in  that  they  measure  the 
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strength  of  the  relationship  between  values  of  the  series 
of  various  lags  after  the  effects  of  other  lags  have  been 
removed.  In  effect/  they  show  the  relative  strength  of  the 
relationship  that  exists  for  varying  time  lags.  For  time 
series  models  of  the  types  to  be  considered,  the  partial 
autocorrelation  coefficients  can  be  calculated  several 
ways.  One  method  of  calculation,  described  by  Box  euid 
Jenkins  [Ref.  4,  Sec.  3.2.5]  begins  with  the  following 
equation,  satisfied  by  the  autocorrelation  function,  where 
the  are  autoregressive  pareuneters: 

Pj  “ ‘‘*klPj-l  + •••  + ♦k^k-lPj-k+l  ♦kkPj-k  ' ^ * 1/2,...,  k 

This  leads  to  what  are  known  as  the  Zule-Walker  equations, 
which  may  be  written 


1 

Pi 

^2  * • * 

Pk-1 

‘"kl 

Pi 

1 

• • • 

Pk-2 

♦k2 

‘>2 

A 

• 

• 
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• 
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• 

• 

• 
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_pk-l 

Pk-2 

Pk-3  •• 

1 

♦ 

kk 

Pk 

Solving  these  equations  for  k » 1,2,3,  ...,  successively , 
we  obt2d.n 


'11 


'1  ' 
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r etc 


In  general,  for  (t)^,  the  determinant  in  the  numerator  has 
the  same  elements  as  that  in  the  denominator,  but  with  the 
last  column  replaced  with  Pj^.  The  quantity  a fiinction 

of  the  lag  k,  is  called  the  partial  autocorrelation  function. 
Other  methods  for  calculating  partial  autocorrelations  include 
successively  fitting  autoregressive  processes  of  orders 
1,  2,  ...,  k by  least  squares  and  picking  off  estimates  of 
^11'  ^22'  ***'  ^kk  lest  coefficient  fitted  at  each 

stage.  Another  method,  discussed  in  i^pendix  A3. 2 of 
Ref.  4 and  due  to  Dtirbin,  generates  estimates  of  an  auto- 
regressive process  of  order  k^-l  recursively  from  estimates 
of  autoregressive  processed  orders  k and  less.  It  is 
derived  by  observing  recursive  relationships  fron  the 
Yule-Walker  equations.  The  recursive  formulas  are: 


j » 1,2,  . . . , k,  and 


4) 


r - 

^k+1  " >i*kj  ^k+l-j 

k+i,k+i  “ ir~ 

" ” ^ *ki 

j-1  ''J  ^ 


Another  way  of  looking  at  partial  autocorrelations 
is  to  consider  a time  series  emd  the  model 

*t  ’ *^0  ^'l^t-l  ^2^2  + •••  + Vl*t-(k-l) 

where  the  b values  are  the  least  squares  estimates  of  the 
linear  regression  coefficients  (6’s)  in  the  model 


*t  ’ ®l*t-l  + •••  + ®k-l*t-(k-l)  ®t  * 

Let  be  the  residual  of  z^  after  removing  the  linear 

effect  of  •••»  *t-k+l  that  “ *t  “ ®t  * 

Then,  the  partial  autocorrelation  of  lag  k,  denoted 

is  defined  to  be  the  simple  autocorrelation  of  lag  k for 

the  adjusted  series  (z^,z-,  z ). 

JL  z n 

The  exact  expression  for  partial  autocorrelations  of 
a moving  average  process  is  complicated,  but  an  approximate 
duality  between  autocorrelations  of  eui  autoregressive 
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process  and  pautos  of  a moving  average  process  in^lies 
that  the  pautos  of  a moving  average  process  would  have  the 
seune  general  behavior  as  the  autocorrelations  of  em  auto- 
regressive process  of  the  S2une  order.  This  approximate 
duality  also  inplies  that  the  pautos  of  an  autoregressive 
process  of  order  p will  have  the  same  general  behavior  as 
the  autocorrelations  of  a moving  average  process  of  order 
p.  ThuSf  one  function  can  be  examined  and  used  to  confirm 
the  other,  emd  may  be  of  use  in  model  identification,  discussed 
in  section  G of  this  chapter. 

G.  COMMENTS  ON  MODEL  SELECTION  EMPLOYING  AUTOCORRELATIONS 

AND  PARTIAL  AUTOCORRELATIONS 

The  purpose  of  this  section  is  to  provide  a brief  over- 
view of  the  methodolo^  for  selecting  the  proper  model  for 
a given  time  series,  through  both  examination  or  prior 
knowledge  of  the  data  itself  and  excunination  of  plots  of 
the  autocorrelation  and  partial  autocorrelation  functions. 

A not  yet  stationary  series  should  first  be  differenced 
or  trauisformed  until  it  appears  to  be  stationary;  character- 
istic of  the  autocorrelation  plot  (correlogram)  of  a non- 
stationary  series  is  a very  slow  damping  out  of  the  auto- 
correlations. When  this  property  of  the  correlogram  is 
observed,  the  emalyst  should  difference  the  series  vmtil  no 
further  improvement  is  reached.  For  ARMA(p,q)  models,  the 
tentative  identification  of  the  model  class,  that  is,  the 
determination  of  p and  q,  can  be  done  by  exaunining  the 


sample  autocorrelation  auid  sample  partial  autocorrelation 
functions  of  the  given  time  series  with  the  theoretical  auto 
and  pauto  functions  of  members  of  the  general  linear  class. 
For  most  stationary  time  series,  an  adequate  fit  can  be 
found  in  a model  with  p and  q relatively  small,  say  three 
or  less. 

The  following  figure  displays  examples  of  autocorrela- 
tion and  partial  autocorrelation  versus  lag  plots  for  several 
basic  model  types,  described  in  Ref.  17. 

Figure  3.  Examples  of  Auto  and  Pauto  Plots 
Figure  3(a).  Auto  and  Pauto  vs.  Lag  for  AR(1)  Model 
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As  cem  be  seen  from  the  figiires,  there  is  a wide  variety 
of  possibilities  for  appearemce  of  the  correlogram  and 
peurtial  autocorrelation  plots  for  even  the  simplest  of 
models.  Sometimes  the  pattern  of  estimated  autos  and  pautos 
can  be  quickly  classified  in  terms  of  one  of  the  sis^ler 
basic  models;  however,  most  real  data  will  generate  auto 
and  pauto  plots  that  resemble  those  of  Figure  3 (e) , those 
of  a mixed  model.  Here  some  free  association  may  be  re- 
quired to  infer  a pattern  from  the  correlogram,  or  more 
than  one  pattern  may  be  implied.  However,  most  of  the  time 
a selection  of  something  is  possible.  As  is  evident,  this 
is  a somewhat  subjective  process,  where  the  quality  of 
identification  will  improve  with  experience. 

The  following  table,  adapted  from  Box  and  Jenkins 
[Ref.  4],  summarizes  the  properties  of  autoregressive, 
moving  average  and  mixed  ARMA  processes;  an  understemding 
of  its  contents  will  be  helpful  as  thumbrules  for  initial 
model  identification. 

It  can  be  noted  from  examination  of  Table  II  that  there 
exists  what  Ccui  be  termed  a "duality"  relationship  between 
the  autocorrelations  and  partial  autocorrelations  of  pure 
autoregressive  and  moving  average  processes.  For  example, 
the  plot  of  autocorrelation  versus  lag  for  an  AR(p)  process 
would  appear  the  same  as  the  plot  of  partial  autocorrela- 
tion versus  lag  for  an  MA(q)  process,  for  p * q.  Similarly, 
the  plot  of  partial  autocorrelation  versus  lag  for  the  AR(p) 
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Table  II.  Properties  of  Basic  ARMA  Processes 

process  would  appear  the  seune  as  the  plot  of  autocorrelation 
versus  lag  for  an  MA(q)  process,  where  p “ q. 

For  seasonal  series,  where  differencing  has  been  finable 
to  remove  all  the  nonstationary  seasonal  behavior,  the  auto- 
correlation euid  partial  autocorrelation  plots  will  generally 
exhibit  a "spike"  at  lags  equal  to  integer  multiples  of  the 
period  of  seasonality.  For  example,  if  the  data  are  seasonal 
with  a period  of  12,  then  the  autocorrelations 
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^12*  ^24'  ^36'  ***  partial  autocorrelations 

*^12  12'  *^24  24'  *^36  36'  ***  ^ amplified.  This  would 

indicate  a need  for  one  or  more  seasonal  p2ureuneters  to  be 
included  in  the  model.  The  following  figure  is  an  illus* 
tration  of  such  an  auto/pauto  plot. 


Figure  4.  Exan^le  Auto  and  Pauto  versus  lag  for  mixed 
ARMA(p,q) , with  seasonality  remaining  of 
period  s » 6. 
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H.  ESTIMATION  OF  MODEL  PAPAMETEBS 

Suppose  a time  series  has  been  tentatively  identified 
as  an  ARIMA(p,d/q)  model: 

(|»(B)  7*^  + 0(B)  a^  . 

There  are  then  p+q+2  unknown  par2uneters  of  the  model  to  be 

2 

estimated,  (4>j^,4»2»  •••»  ®o' ®1' ®2'  ***'  ®q'^a^  * 
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Box-Jenklns  procedure  separates  this  estimation  into  two 

parts.  First,  estimates  are  obtained  for  the  AR  and  MA 

parameters  <() » •••»  <l>p)  9 * ***'  ' 

2 

then  estimates  are  calculated  for  9 auid  o , which  are  func- 

O CL 

tions  of  <|>  and  9.  . 

The  usual  procedure  is  to  select  those  parameter  values 
$ and  9 that  minimize  the  sum  of  squared  model  errors 
(residuals).  Let  “ *^*2  ~ ***  ” *^p^ 

W.  - V^(z.  - u ,)  • Then,  it  can  be  shown  that  u « EtW  ] 
and  that  the  model  can  be  rewritten  as 

(W^  - u^)  * 9(B)  a^  , 


or 


= 9“^(B)  (f(B)  (W^  - u^)  . 


Now,  let  W,  the  sample  mean,  be  the  estimate  of  (if 
d > 0,  then  y^  is  usually  zero) . Also,  let 


a^  = 9"^(B)  $(B)  (W^  - W) . 


\ 

: 

5 

•j 


Set 

N 

s(i,i)  = I ' 

t«l 

where  N - n-d  is  the  length  of  the  differenced  series.  The 
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objective  is  to  select  those  pareuneters  £ and  £ such  that 

S(^,9^)  is  minimized  over  all  values  of  the  (p-H}) -dimensioned 

pareuneter  space.  Since  the  equation  in  S is  nonlinear 

in  the  parameters/  iterative  search  methods  are  generally 

used  to  determine  optimal  model  parameters.  Having  found 

2 

these  optimal  parameters  and  the  minimum  S,  a C2tn  be 

cl 

estimated  using 

52  . s(|/§)/N-p-q  . 


Chapter  III  contains  a more  detailed  description  of  a 
nonlinear  algorithm  to  estimate  pareuneters  of  a generalized 
seasonal  or  nonseasonal  ARIMA  model. 

I.  DIAGNOSTIC  CHECKING  OF  MODELS 

After  the  model  has  been  tentatively  identified  and 
par2uneter  estimates  have  been  calculated,  the  next  task  is 
to  test  whether  or  not  the  original  model  specification 
was  correct  and  that  the  model  itself  is  adequate  in  fore- 
casting power.  The  process  of  testing  the  model  cam  teUce 
many  forms,  but  will  usually  include  at  least  the  following 
two  steps : 

1.  Generate  a simulated  series  from  the  estimated 
model  and  compare  the  simulated  series  and  its 
auto  and  pauto  functions  with  the  original  series 
and  its  respective  auto  and  pauto  functions.  This 
comparison  is  essentially  subjective. 
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2.  Calculate  the  residuals  of  the  estimated  model,  the 
a^'s,  and  compare  the  properties  of  the  residuals 
with  those  assumed  for  the  shock  terms  of  the  actual 
process.  The  residuals  should  be  normally  distri- 
buted and  uncorrelated  with  each  other;  that  is,  there 
should  be  no  discernable  structure  in  the  residuals. 
There  are  m2my  quantitative  statistical  tests  and 
data  analytic  tools  that  cam  be  applied  to  the 
residuals  to  test  hypotheses  of  normality  and  zero 
autocorrelation . 

A plot  of  the  autocorrelation  and  partial  autocorrelation 
functions  of  the  residuals  will  provide  not  only  a test  of 
whether  or  not  the  residuals  are  uncorrelated,  but,  if  they 
are  correlated,  the  plots  can  be  used  to  suggest  improve- 
ments to  the  model.  For  example,  suppose  the  model  was 
tentatively  specified  as  the  ARMA(1,1)  model  below: 

(1  - 0.5B) (z^  - 2)  » (1  + 0.7B)  a^  . 

Suppose  also  that  the  autos  and  pautos  of  the  model  residuals 

suggested  the  model  (1  - 0.3B)  a^  = u^,  where  the  u^'s 

2 

are  white  noise,  uncorrelated  with  variance  a^.  Then, 
these  two  models  can  be  combined  into  the  ABMA(2,1)  model: 

Cl  - 0.3B)  (1  - 0.5B)  (z^  - 2)  » Cl  + 0.7B)  u^  . 

This  should  be  a refinement  over  the  original  model. 
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J.  FORECASTING 

The  objective  in  forecasting  is  to  predict  future 
values  of  the  time  series  with  as  little  error  as  possible. 
The  criterion  most  often  used  for  forecasting  is  to  calcu- 
late that  forecast  which  minimizes  the  expected  meem  square 
forecast  error.  Therefore,  if  z^(i)  denotes  the  forecast 
for  lead  time  I from  time  origin  t of  the  value 
the  objective  is  to  find  z^(i)  such  that  the  objective 
function 


E[(Zt+j^  - z^(i)  ] 

is  a minimum.  This  forecast  is  given  by  taking  z^d)  as 
the  conditional  expectation  of  given  z^,  z^; 


The  forecast  can  be  easily  generated  recursively  from  the 
mathematical  model,  utilizing  the  facts  that 


E [z^_ j 1 z^, . . . , Zil  “ *t— j 


for  j - 0,1,2...  (where  t is  the  current  time),  emd 


E[Zt+j |Zt/ • • • /^il  * for  j « 1,2, 


and 


0 for  times  > t 

E[a^]  - 

a^  for  times  ^ t 

For  example,  consider  the  Box** Jenkins  model 
(1  - 0.5B  + 0.6B^)  2^  - (1  + 0.3B) 

which  could  be  written  in  the  exp2mded  form  as 


0.52^_j^  - 0.62^_2 


‘t-1 


Now,  assume  it  is  known  that  z^qq  is  1.4,  Zg^  is  1.0,  and 
the  calculation  has  been  made  to  obtain  the  residual 
a^QQ  » 0.2.  Then  the  forecasts  of  ^102  ^103 

made  from  an  origin  of  t = 100  can  be  calculated  from  the 
model  as  follows: 

first,  for  4*1,  set 


*100^^^  “ ®^®10ll ^100'*99' * * * '^1^  ' 


^100^^^ 


EtO.SZioo  “ *101  °*^*1Q0^ 


and 


ZiqqH)  - (0.5)  (1.4)  - (0.6)  (1.0)  + 0+  (0.3)  (0.2)  ; 


then  our  forecast 


. 0.7  - 0.6  + 0.06  » 0.16  ; 


Next#  for  I • 2,  set 


*100^^^  “ ®^*102l*101'*100'*99' • '*1^  ' 


*100^^^  * " °*®*100  ^102 


*100'^99'***'*1^  ' 


£^Oo(2)  = (0.5)(S^0q(1))  - (0.6)(z^00)  + 0 + 0 , 


ZlOo(2)  = (0.5)  (0.16)  - (0.6)  (1.4)  - -0.76  ; 


similarly,  for  A * 3,  set 


*100^^^  “ ®^*103l*102'*101'*100,*99"* * * '*1^  ' 


*100^^^  * ®t°*5Zio2  “ °*®*101  *103  °*^*102 


*100'*99'--*'*1^ 


SioqO)  - (0.5)  (S^qq(2))  - (0.6)(z^qq(1))  + 0 + 0 , 


/ 
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emd  then 


^ ; 

! (* 


: t 


i 

I 5" 


! I I 


(0.5) (-0.76)  - (0.6) (0.16) 


-0.48  . 


Cleeurly/  this  process  can  be  continued  into  the  future  as 
long  as  desired,  recognizing  that  the  expected  forecast 
accuracy  will  decrease  as  I increases. 

Let 


e^(i) 


(zt^,  - 


denote  the  forecast  error  for  lead  time  I beyond  the  forecast 
time  origin  t.  It  cam  be  shown  that  e^(ll)  is  given  by 
the  relation 


et(i)  - '^i-lS+1  ' 


where  the  scalar  weights  are  determined  from  the  equation 


i|»(B)  - r^(B)  (1  - 8)“*^  0(B) 


using  the  known  parauneter  values.  The  variamce  of  the 
forecast  error  is  given  by 


i-l 


E((eJ()l)]  - (1  + il*!  + 


{1  + I 


j-i 


'j^'^a 
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From  these  relationships,  a confidence  interval  of  n 
standard  deviations  about  a forecast  of  lead  time  i would 
be  given  by: 


Cn(2t(4)) 


± n[l  + 


I '•'jl 


j-1 


The  value  of  n could  also  be  though  of  as  the  percentage 
point  for  the  desired  l^el  of  confidence  using  the  Normal 
distribution  [Ref.  4].  It  can  be  noted  from  the  definition 
of  s^(£.)  edsove  that  the  one-step-ahead  forecast  error, 
e^(l)  is  simply  the  value 


®t*i’  * “t+i 


This  is  an  explanation  for  the  use  of  the  term  "residual" 
to  refer  to  the  white  noise  or  shock  terms.  Additionally, 
from  the  foregoing  it  is  evident  that  the  forecast  error 
variance  is  a non-decreasing  function  of  the  length  of 
the  forecast  lead  time  i;  therefore,  the  confidence  bands 
must  become  wider  as  the  forecast  lead  time  increases. 


50 


III.  GENERALIZED  ARIMA  MODEL  PARAMETER  ESTIMATION 


The  purpose  of  this  chapter  is  to  describe  the  compu- 
tational details  of  the  algorithm  WMARQRDT,  employed  with 
appropriate  subroutines  to  estimate  Box-JenJcins  pareuneters 
of  a generalized  seasonal  (or  non-seasonal)  ARIMA  model. 
This  generalized  model  can  be  written  parsimoniously  as 
follows: 


<|»(B)*(B®)W^  - 0Q+0(B)0(B®)a^  , 


where 


"t 

emd 


z^  is  the  tramsformed  time  series  (for  example, 
z^  * In  z^  or  z^  « /z^)  . 

The  model  employs  a Marquafdt-type  (see  Ref.  13]  non- 
linear least  squaures  algorithm  for  determination  of  the 
model  parameters,  searching  in  the  parauneter  space  for  the 
set  of  paurameters  which  minimizes  the  sum  of  the  squared  model 


residuals. 


A.  MATHEMATICAL  DISCUSSION  OF  THE  MARQUARDT  ALGORITHM 

Reference  13  describes  Marquardt's  methodology  for 
development  of  a general  algorithm  for  least-squares  esti- 
mation of  nonlinear  parameters.  Most  algorithms  for  least- 
squeures  estimation  of  nonlinear  parameters  have  been  cen- 
tered 2d}Out  either  a pure  linear  iterative  model  based  on 
a Taylor  series  expansion  or  using  s<xne  form  of  the  method 
of  steepest  descent  (maximum  negative  gradient) . Since 
both  methods  have  severe  potential  pitfalls  (Taylor  series 
due  to  divergence  of  the  iterates  due  to  an  unfortunate 
choice  of  initial  values/  the  steepest  descent  due  to  slow 
convergence  after  the  first  few  iterations) , a method  termed 
"Maximum  neighborhood"  method  was  developed  by  Marquardt 
which,  in  effect,  performs  an  optimum  interpolation  between 
the  Taylor  series  method  and  the  method  of  steepest  descent, 
the  interpolation  being  based  upon  the  meocimum  neighborhood 
in  which  the  truncated  Taylor  series  gives  an  adequate 
representation  of  the  nonlinear  model. 

As  discussed  in  Ref.  14,  the  idea  of  Marquardt's  method 
can  be  briefly  explained  as  follows.  Assume  the  algorithm 
begins  with  a vector  of  initial  parameters  * 

If  the  method  of  steepest  descent  is  applied,  a certain 
vector  direction  h^,  where  the  subscript  g represents  the 
gradient  (which  gives  the  direction  in  which  the  rate  of 
cheuige  is  the  greatest) , is  obtained  for  movement  away  from 
the  initial  point  values  in  the  parameter  space.  Due  to 
potential  nonlinearities  in  the  sum-of-squared-residuals 
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surface,  S 0^, 0^) , in  the  parameter  space,  may 

be  the  best  local  direction  in  which  to  move  toward  opti- 
mality (minimal  sum  of  squared  model  residuals) , but  may 
not  be  the  best  direction.  However,  the  best  direction 
must  be  within  90  degrees  [Ref.  8]  of  h , or  else  the 
value  of  S (4) , 9 , 0)  will  increase  locally.  Now,  the 
linearization,  or  Taylor  series  method,  may  result  in  a 
different  correction  vector,  h^.  Marquardt  foxind  experi- 
mentally that  for  a number  of  practical  problems,  the  angle, 
say  C,  between  the  vectors  h^  and  h^,  fell  in  the  range 
80  degrees  < C < 90  degrees.  In  other  words,  the  two 
directions  were  nearly  normal  in  most  cases.  The  Marquardt 
algorithm  then,  provides  a method  for  interpolation  between 
the  vectors  h^  and  h^,  employing  a sutiable  step  size,  which, 
in  general,  converges  more  quickly  than  either  of  the  two 
"parent"  methods  alone.  The  remainder  of  this  section  des- 
cribes the  mechanics  of  the  algorithm,  as  applied  to  the 
estimation  of  Box-Jenkins  parameters. 


B.  MODEL  INPUT  INFORMATION 

The  following  is  the  required  input  information  to 
the  WMARQRDT  parameter  estimation  model; 


the  transformed  time  series  (for  example, 
= In  z^,  or  * /z^ 

the  length  of  the  time  series 

the  length  of  the  seasonal  period  (for  no 
seasonality,  s =*  1) 


Prior  to  initiating  the  sum  of  squares  calculation  pro- 
cess, the  model  unravels  the  polynomial  of  AR  and  MA  operators 


and  forms  a vector  of  coefficients  for  both  the  AR  and  MA 
sides  of  the  model*  For  example,  consider  the  relatively 
simple  (2,1)  x (2,1) ^2  ARIMA  model,  written 

(1-<D^B-(|)2B^)  (1  - • (1-6^8-028^)  (l-e3^  jL2®^^^“t* 

Unraveling,  the  model  becomes  first: 

and  then  takes  the  form: 

\“*1, 12\-12"‘^l\-l'^‘*‘l*l,  12\-13'*2\-2‘^*2*1 , 12’'t-14 

“ *t"®l,12S-12"®lS-l'‘’®l®l,12S-13"®2^t-2'*’®2®l,12S-14  . 

This  is  represented  in  the  computer  as  four  vectors:'  two 
vectors,  tt  and  y,  composed  of  ^e  expemsion  AR  emd  MA 
coefficients  respectively,  and  two  vectors  composed  of  the 
indices  of  the  time  series  terms  associated  with  the  coeffi- 
cients. For  the  model  in  the  example  eUsove,  these  vectors 
would  be  written  as  follows: 

TT  « ^^'*1,12'***!' ^'**1^^'*1,12^ '***2' ^*^2^  ^*1,12^  ^ 

INDEX^  * (0,  12,  1,  13,  2,  14) 

Y » ^^'®1,12'®1' ^®1^  ^®1,12^ '®2' ^®2^^^,12^  ^ 


i: 


This  format  for  iinraveling  of  a generalized  product  of 
polynomials  of  seasonal  and  non-seasonal  pareuneters  permits 
simplified  calculations  in  the  determination  of  the  sum  of 
squared  residuals,  since  the  eurithmetio  need  only  be  per- 
formed for  the  non-zero  parameter  values,  those  contained 
in  the  it  and  y vectors.  In  the  actual  calculations,  per- 
formed in  FORTRAN,  the  index  vector  origin  is  translated 
to  the  right  by  the  amount  p + Ps  + q -«•  Qs,  called  lAODIT 
in  the  model,  in  order  to  prevent  obtaining  zero  or  negative 
subscripts . 

The  calculation  of  the  residual  sum  of  squares  is  a 
three  stage  process.  First,  the  model  calculates  ten 
values  of  e^,  the  "forward"  white  noise  terms,  using  the 
following  relationship: 

p+Ps  q+Os 

•t  * "t  • ^ 'iVi  * ^ 

i»2  j-2 

The  assumption  is  made  that  equal 

to  zero.  Having  found  the  desired  values  for  (e^Q,eg, . . . ,e^) , 
we  C2m  now  backcast  the  necessary  values  of  W^,  in  order 
to  enable  us  to  calculate  the  required  white  noise  terms 
a^,  as  follows: 
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p+Ps 

*t  * ^ ’iVt 


q+Qs 

j-2 


solving  for 

W^'s  backcast,  we  can  proceed  to  calculate  the  estimated 
white  noise  terms  (residuals)  using  this  relationship: 


p+Ps  q+Qs 

I + I Yj^t-j 

i-2  j-2 


The  values  of  a^  are  calculated  for  t » 1,  ...,  N.  Once  the 
a^’s  are  calculated,  it  is  a simple  matter  to  calculate  the 
stun  of  squares,  using 


S (<^,  $,9,0) 


V ^ 
^ ^t 


Having  calculated  the  residual  sum  of  squares  for  the 
initial  parameter  estimates,  we  form  the  parameters  into 
a vector  0,  such  that 


(0^,02^  •••/  0j^)  t 


where  k ■ p+P+q+Q  ; that  is. 


®q  s ' * * ^Q,  s^ 


ovttr  *11  valuAS  of  t * and  1 « Using  the 

residuals  calculated  earlier / the  derivatives  are 

estiaated  numerically  using  a perturbation  of  each  parameter 
of  amount  5,  such  that 


^ • {a^  ( / . • . f Sj^)  “a^  ( ^ f S j^+6 » . . . f Bj^)  J/S  > 

Then  with  {a^}  amd  ^ supplied  for  the  current  parauneter 
values,  the  following  quamtities  are  formed: 


1.  The  kxk  matrix  A = • where 


t®l 


2.  the  vector  G * {g^}  / where 


t»l 


3.  the  scaling  quauitities  = /A^ 


Then,  the  modified  (scaled  and  constrained)  linearized 


correlation  matrix  for  the  model  pareuneters.  Those  equa 


tions  are  then  solved  for  h*,  which  is  scaled  back  to  give 


Then,  the  new  parameter  values  are  formed  from  the  h vector. 


such  that 


and  the  sun  of  squared  residuals  for  the  set  of 

parameters  is  calculated. 

Now,  if  < S(6qjjj),  the  pareuneter  corrections  h 

are  tested.  If  all  elements  of  the  h vector  are  smaller 
than  some  epsilon,  say  .0001,  then  convergence  is  assumed 
and  the  kxk  matrix  is  used  to  calculate  the  covari- 
ance matrix  of  the  estimates;  otherwise,  is  reset  to 

^NEW'  constemt  ^ is  reduced  by  a factor  F,  say  .1,  emd 
computation  returns  to  the  calculation  of  a new  set  of 
derivatives. 

However,  if  the  constraint  parameter 

^ is  increased  by  the  factor  F,  and  confutation  is 
returned  to  the  stage  where  the  A matrix  is  formed.  In 
all  but  xinusual  cases,  a minimal  sum  of  squares  will  be 
found.  However,  an  upper  bound  can  be  placed  on  and 
when  it  is  exceeded,  the  search  is  terminated.  When  con- 
vergence is  reached,  the  desired  output  information  is 
calculated. 


D.  MODEL  OUTPUT  AND  DIAGNOSTIC  CHECKS 


The  output  includes: 


1.  B “ (B^»»«»»Bj^) 


the  least-squares  estimates  of 


<<> 

0 

0 


1 

1 

1 

1 


F 


F 
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2.  a! 


3.  V 


4.  Sj 


5.  R 


ij 


6.  0. 


the  residual  variauice,  ^p,gj^p_Q  S($,$,§,0) 


the  covariance  matrix  of  the  parameter 
estimates,  formed  from 


V - {V^j}  - (A^A)"^ 


the  standard  errors  of  the  peurameter 
estimates , 


Sj^  « /V^  , i » 1, . . . ,p+q+P+Q 


the  correlation  matrix,  obtained  from 


the  overall  constant  term,  where 


9^  = y G , for 
o 


P P . 

(1  - I (1  - I and 

i-1  j=l 


N 

t-1 


I 


Finally,  diagnostic  checks  on  the  model  are  performed. 

This  includes  calculation  of  residual  autocorrelations 
I obtained  from  these  fomulas: 

i I 


rjj  (k) 


cjj(k)/cjj(0) 
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IV.  DESCRIPTION  OP  THE  TIME  SERIES  EDITOR 


In  this  chapter,  descriptions  are  given  for  each  of 
those  progr^uns  in  the  Time  Series  Editor  which  interact 
with  the  user.  There  are  fifteen  separate  program  modules 
currently  included  in  the  Editor;  there  are  other  programs 
utilized  by  the  Editor  that  are  completely  transparent 
to  the  user,  and  are  not  described  here.  However,  listings 
for  all  programs  contained  in  the  Editor  are  provided  in 
Appendix  D to  this  report. 

No  attempt  is  made  in  this  chapter  to  describe  the 
actual  mathematical  calculations  or  algorithms  that  are 
performed  by  the  progrcuns.  Rather,  the  objective  is  to 
provide  general  descriptions  of  what  each  program  can  do 
for  the  user,  and  how  the  user  interacts  with  the  programs. 
The  programs  described  in  this  chapter  are,  in  the  order 
presented:  TIMESER  EXEC,  ZFORMAT,  CMSWORK,  TRANS,  DIFF, 

PLOT,  AUTO,  ESTIMATE,  YESTSEAS,  WMARQRDT,  XSUMSQ,  FORECAST, 
ROOTS,  SIMULATE  and  GENERATE,  and  HELP. 

A.  THE  EXECUTIVE  PROGRAM 

The  heart  of  the  Time  Series  Editor  is  a master  program 
called  TIMESER  EXEC  which  provides  file  control  for  all  of 
the  other  progr2un  modules,  controls  input  amd  output,  super- 
yises  the  necessary  CP/CMS  protocol  and  provides  instructions 
interactively  to  the  user  concerning  the  contents  of  the 


Time  Series  Editor  and  how  each  program  can  be  used. 

The  TIMESER  EXEC  program  is  written  in  the  IBM-360  CP/CMS 
Executive  language.  It  is  the  only  program  in  the  package 
not  written  in  FORTRAN;  consequently,  it  should  be  the  only 
program  that  would  need  modification  if  the  Editor  were  to 
be  adapted  to  euiother  FORTRAN-capable  time-sharing  system. 

After  the  user  has  logged  into  CP/CMS  and  linked  to 
the  disc  space  containing  the  Time  Series  Editor  (instruc- 
tions for  this  procedure  are  provided  in  detail  in  the 
User's  Guide  to  the  Time  Series  Editor,  included  as  Appendix 
A to  this  report) , the  entire  Editor  package  is  made  avail- 
eUsle  to  the  user  by  the  command  TIMESER.  On  entry  of  this 
comand,  the  executive  routine,  TIMESER  EXEC,  can  provide 
a guided  tour  through  the  Editor.  It  will  briefly  des- 
cribe what  the  Editor  can  provide  emd  asks  the  user  what 
option  he  wants  to  use.  On  the  basis  of  the  user's  response, 
the  Editor  then  advises  the  user  what  input  data  is  required 
and  how  it  should  be  entered,  either  through  a data  file 
entered  offline  or  parametric  data  entered  interactively 
at  the  terminal.  When  the  user  selects  an  option  for  exe- 
cution, the  TIMESER  EXEC  routine  loads  the  appropriate  pro- 
gram package  emd  automatically  manipulates  the  required 
input  and  output  data  files.  For  the  more  experienced  user 
who  does  not  require  detailed  user  instructions,  there  is 
a shortened  version  of  the  TIMESER  EXEC.  Upon  entry  into 
the  shortened  version,  the  user  is  immediately  asked  which 
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option  he  desires.  This  version  of  the  Editor  is  entered 
by  adding  the  argument  "s"  when  logging  into  the  Editor; 
that  is,  type  TIMESER  S when  logging  in.  The  TIMESER  S 
(shortened)  version  of  the  Editor  provides  exactly  the 
same  program  options  to  the  user  as  the  longer  TIMESER 
version. 

B.  DATA  AND  FILE  MANAGEMENT 

Whenever  data  is  required,  the  user  is  prompted  at 
the  terminal  by  either  the  TIMESER  EXEC  or  the  program 
module  being  executed.  In  most  cases,  the  necessaury  user 
response  is  a short  alphanumeric  character  input  during 
execution  using  the  terminal  keyboard.  However,  when  the 
time  series  data  itself  is  entered,  the  user  may  enter  the 
data  by  card  deck  offline.  Similarly,  most  output  is  pro- 
vided right  at  the  user's  terminal.  Output  such  as  listings 
of  transformed  series  or  plots  are  sent  to  the  offline 
printer  to  conserve  time  and  provide  hard  copies  of  the 
results.  Detailed  descriptions  of  the  input  euid  output 
requirements  of  each  program  module  are  provided  in  this 
chapter,  and  also  in  abbreviated  tabular  form  in  Appendix 
A,  the  User's  Guide.  However,  there  are  some  general  prin- 
ciples of  data  mzmagement  that  apply  to  all  programs. 

These  are  described  in  this  section. 

1.  Offline  Data  Entry 

The  Time  Series  Editor  requires  that  the  user's 
time  reside  in  FILE  FT02F001.  The  data  can  be  read  offline 
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via  cards  and  directed  to  FILE  FT02F001.  An  example  of  the 
required  format  for  the  input  card  deck  for  the  Time  Series 
Editor  is  provided  in  Figure  5.  It  is  noted  here  that  the 


data  format  expected  by  all  programs  is  FORMAT(5F15.6) , 
preceded  by  a card  containing  the  length  of  the  time  series 
in  FORMAT (13) . 

2.  ZFORMAT  Program 

In  many  cases,  the  user  will  want  to  analyze  time 

series  data  that  are  already  available  ina  format  other 

than  that  required  by  the  Editor.  To  spare  the  user  the 

tedious  task  of  retyping  his  deck  in  the  required  format, 

the  ZFORMAT  program  was  written.  This  program  converts  a 

time  series  data  file  (FILE  FT03001)  written  in  any  FORTRAN 

format  into  the  proper  FORMAT (5F15. 6)  onto  FILE  FT02F001. 

/ 

The  length  of  the  series  is  entered  as  an  input  on  the 
terminal.  The  data  are  then  ready  to  use  in  the  Editor, 
with  the  original  file  left  in  FILE  FT03F001.  Figure  6 
shows  a seunple  input  deck  for  such  non-standard  data. 

Figure  7 gives  a saunple  session  where  ZFORMAT  is  used. 


igure  5.  Card  deck  arrangement  for  OFFLINE 
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timesec  s 

ENTER  LETTER  FOR  OPTION  YOU  WANT, 

z 

IS  YOUR  DATA  IN  FILE  FT03F001? 

y 

EXECUTION  BEGINS... 

IS  YOUR  TIME  SERIES  DATA  NOW  IN  FILE  FT03F001? 


ENTER  THE  LENGTH  OF  YOUR  TIME  SERIES  VIA  FORMAT  13. 

010 

NOW  ENTER  THE  FORTRAN  FORMAT  FOR  YOUR  TIME  SERIES  DATA, 
INCLUDING  PARENTHESES;  FOR  EXAMPLE,  TYPE:  (5F15.7). 
flO.2 

YOUR  FORMAT  IS: 

FORMATF10.2 
IS  THIS  CORRECT? 


NOW  ENTER  THE  FORTRAN  FORMAT  FOR  YOUR  TIME  SERIES  DATA, 
INCLUDING  PARENTHESES;  FOR  EXAMPLE,  TYPE:  (5F15.7). 
(flO.2) 

YOUR  FORMAT  IS: 

FORMAT(F10.2) 

IS  THIS  CORRECT? 


YOUR  DATA  IS  NOW  TRANSFORMED  INTO  THE  PROPER  FORMAT 
FOR  USE  IN  THE  TIME  SERIES  EDITOR, LOCATED  IN  FILE 
FT02F001.  YOUR  ORIGINAL  DATA  IS  IN  FILE  FT03F001. 

DO  YOU  WANT  TO  GO  AGAIN? 
n 

CONTROL  RETURNED  TO  CMS 

R; 


Figure  7.  Sample  user  session  with  ZFORMAT  program 
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3.  Interactive  Keyboard  Input 
a.  Numeric  Input 

When  the  user  is  prompted  at  the  terminal  to 
enter  such  numerical  values  as  the  number  of  AS  par^uneters, 
estimates  of  pareuneters  or  the  number  of  differences  to  be 
taken,  he  must  enter  these  values  before  program  execution 
can  continue.  The  user  should  enter  the  data  according  to 
stemdard  FORTRAN  practice.  For  example,  integer  data  should 
be  entered  (without  a decimal  point)  for  counts,  such  as 
the  length  of  a seasonal  period,  and  for  names  beginning  , 

with  the  letters  I through  N?  floating-point  data,  such  as 
initial  estimates  of  parameters,  should  be  entered  with  a 
decimal  point.  Because  a typed  decimal  point  overrides  a 
floating-point  format,  it  is  not  necessary  for  the  user 
to  concern  himself  with  the  format  for  floating-point  input 
data.  However,  care  must  be  taUcen  when  entering  integer 
data  because  it  must  be  right-justified  in  its  format  field. 

The  user  is  advised  in  each  case  when  the  format  is  other 
than  II.  For  example,  suppose  the  user  desires  to  perform 
analysis  on  a series  having  only  90  observations.  If  the 
program  that  he  is  executing  requires  the  length  of  the  time 
series,  the  user  will  receive  the  following  request: 

ENTER  LENGTH  OF  THE  TIME  SERIES  VIA  FORMAT  13 


The  user  would  then  enter: 
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colunm  123 


Cwhere  b represents  a blank  space) 


If  the  blank  were  omitted,  the  program  would  read  the 
series  length  as  900,  and  problems  in  execution  would 
occur. 

b.  AlpheUsetic  Input 

The  Editor  often  asks  the  user  to  respond  with 
alphabetic  input.  Generally,  this  is  in  the  form  of  the 
response  to  a question  requiring  a yes  or  no  answer,  or  a 
response  to  the  question  of  which  program  to  execute.  The 
editor  has  been  programmed  to  read  only  the  first  letter 
of  such  responses.  Therefore,  the  user  need  only  enter  the 
first  letter  of  the  response  for  each  such  inquiry.  For 
example,  the  user  should  enter  Y for  yes,  N for  no,  P for 
the  PLOT  option,  X for  the  XSUMSQ  option,  W for  the 
WMAPQRDT  option,  etc.  The  other  alpheibetic  input  is 
generally  entry  of  titles  to  plots;  these  can  be  any 
combination  of  numeric  euid  alphedsetic  characters,  as  long 
as  they  do  not  exceed  72  columns  in  length.  For  example,  a 
suitable  title  would  be: 


PLOT  OF  AUTOS  AND  PAUTOS  FOR  SERIES  C DATA  (2  ORDINARY  DIFFERENCES) 


4.  Output 

Most  of  the  results  are  written  out  at  the  user^s 
terminal.  However,  in  some  cases  the  output  is  printed 
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offline  to  conserve  time.  Such  results  as  plots  and 
transformed  time  series  are  written  onto  various  files 
(PT03F001,  FT09F001,  FT08F001) , and  the  TIMESER  EXEC  pro- 
gram causes  them  to  be  automatically  printed  offline  under 
the  user's  Identification  number  (USERID).  The  plots  and 
files  can  also  be  printed  out  at  the  terminal  using  the 
usual  CMS  commands  In  the  CMS  mode,  outside  the  Editor,  or 
In  the  CMSHORK  program  (see  next  section) . However,  the 
user  should  be  aware  that  plotting  graphs  at  the  terminal 
Is  generally  quite  a slow  process. 

5.  CMSWORK  Program 

This  program  permits  the  user  to  enter  the  CMS 
environment  to  perform  routine  CMS  functions  without  leaving 
the  Time  Series  Editor.  Commonly  used  CMS  commands  Include 
file  neune  alteration,  file  erasure,  file  offline  printing 
or  punching,  obtaining  a disc  status,  and  listing  all  disc- 
resident  files.  However,  essentially  all  CMS  commands  ceui 
be  Issued  while  in  CMSWORK.  Figure  8 provides  an  example 
of  the  use  of  the  CMSWORK  program. 

C . TRANS  PROGRAM 

The  TRANS  program  takes  a given  time  series  In  File 
F102F001  and  performs  a transformation  upon  each  data  ele- 
ment. The  options  Include  a change  of  scale,  a power  trans- 
formation, a shift  of  the  origin,  and  a natural  logarithm 
treuisformation.  The  original  data  is  written  onto  a file 
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ENTER  LETTER  FOR  OPTION  YOU  WANT. 

C 

ENTER  DESIRED  CP/CMS  COMMANDS,  ONE  PER  LINE. 

WHEN  FINISHED  TYPE:  iGOTO  -QUES 
stat 

P (191):  30  FILES;  263  REC  IN  USE,  33  LEFT  (of  296),  89%  FULL  (2  CYL) 
listf  * ft02f002 
FILE  NOT  FOUND 
17.14.13  LISTF  * FT02F002 
tl!  E(00002)  Ml 
listf  * ft02f001 

FILENAME  FILETYPE  MODE  NO. REC.  DATE 


SERG 

FT02F001 

PI 

3 

9/16 

SERC 

FT02F001 

PI 

5 

9/16 

LNSERG 

FT02F001 

PI 

3 

9/16 

FILE 

FT02F001 

PI 

3 

9/16 

cp  q f 
FILES:- 

NO  RDR,  NO 

PRT, 

NO  PUN 

erase  file  ftOSfOOl 
offline  print  file  ft02f001 
offline  print  serg  ft02f001 
offline  print  sere  ft32f001 
erase  file  ft02f001 
stat 

P (191):  20  FILES;  240  REC  IN  USE,  56  LEFT  (of  296),  81%  FULL  (2  CYL) 
&goto  -ques 

DO  YOU  WANT  TO  GO  AGAIN? 
n 

CONTROL  RETURNED  TO  CMS 

R; 


Figure  8.  Sanple  user  session  with  CMSWORK  program 
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called  DATA  FT02F001,  while  the  treuisfomed  series  is 
written  onto  PILE  FT02F001.  The  peurzuneters  of  the  trans- 
form itself  are  written  onto  FILE  FT07F001.  The  TRANS 
program  is  self-contained,  using  no  other  programs. 


D.  DIFF  PROGRAM 

The  DIFF  program  teUces  a given  time  series  in  FILE 
FT02F001,  and  performs  non-seasonal  and/or  seasonal  differ^ 
encing  operations  on  the  series  in  order  to  help  achieve 
station arity  for  proper  modeling.  The  user  must  indicate 
if  seasonal  differences  are  required;  if  so,  the  length  of 
the  season  is  input.  The  user  is  then  required  to  enter 
the  number  of  seasonal  auid  nonseasonal  differences.  After 
execution  of  DIFF,  the  original  data  are  in  FILE  FT02P001, 
and  the  differenced  data  are  written  onto  FILE  FT03F001. 

To  use  the  differenced  data  in  any  other  progr^,  CMSWORK 
prograua  must  be  entered  and  the  files  altered  to  save  the 
original  data  (if  desired)  and  to  transfer  the  differenced 
data  to  FILE  FT02F001.  The  CMS  commands  to  perform  this 
shift  of  files  would  look  like: 


ALTER  FILE  FT02F001  Pi  SAVEDATA  FT02F001  Pi 
R; 

ALTER  FILE  FT03F001  Pi  PILE  FT02F001  Pi 
R; 


The  DIFF  program  uses  the  IMSL  s^lbroutine  called  FTRDIF 
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E.  PLOT  PROGRAM  ^ 

The  PLOT  program  plots  any  given  time  series  which 
resides  in  FILE  FT02F001.  Other  them  the  time  series, 
which  must  be  on  the  user's  disc  or  entered  offline,  the 
program  requires  only  that  an  identification  title  for  the 
plot  itself  be  entered  during  execution.  The  plot  is  auto- 
matically printed  offline.  The  PLOT  progreun  uses  a modified 
PLOTP  subroutine  called  PLOTS  from  the  IBM  Scientific 
Subroutine  Package  Library  (SSPLIB) . 


F . AUTO  PROGRAM 

The  AUTO  program  takes  a given  time  series  in  FILE 
FT02F001  emd  calculates  summary  statistics  of  utility  in 
Box- Jenkins  modeling.  The  statistics  include  the  sample 
mean,  variance,  autocorrelations  and  peurtial  autocorrela- 
tions for  lags  one  though  40.  Part  of  those  statistics 
are  printed  out -at  the  user's  terminal,  while  plots  of  the 
autocorrelations  and  partial  autocorrelations  versus  lag 
are  printed  offline.  These  plots  will  also  be  written  onto 
the  user's  disc  in  FILE  FT08F001,  amd  may  be  plotted  at 
the  teinninal  in  CMS  outside  the  Editor  or  in  the  CMSHORK 
program;  recall  that  this  is  a time-consuming  process.  The 
AUTO  program  provides  the  user  with  essential  information 
cUsout  stationarity,  seasonality  and  model  identification. 
With  the  Box- Jenkins  procedures,  the  second  moments  (auto- 
correlations emd  partial  autocorrelations)  are  valuable 
tools  in  tentative  model  identification.  Additionally,  AUTO 
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is  useful  in  the  diagnostic  checkout  phase  of  model 
building*  when  examining  the  model  residuals  for  structure. 

The  AUTO  progreun  uses  the  IMSL  subroutine  FTAUTO. 

G.  ESTIMATE  PROGKAM 

After  the  user  has  tentatively  Identified  a model  for 
a series  through  euialysis  of  the  series  plot,  autocorrela- 
tions auid  partial  autocorrelations,  ESTIMATE  should  be 
executed  to  calculate  maximum  likelihood  estimates  of  the 
model  parameters.  ESTIMATE  is  written  to  calculate  parameters 
for  non-seasonal  models  only;  for  seasonal  models,  the 
programs  YESTSEAS  and  WMARQRDT  should  be  used  for  parameter 
estimation.  The  ESTIMATE  program  requires  that  the  series 
to  be  modeled  (already  tremsforroed  and  differenced  to 
achieve  stationarity)  reside  in  FILE  FT02F001;  additionally, 
the  prograun  will  instruct  the  user  to  input  the  number  of 

i 

MA  and  AR  parcuneters  in  the  model,  as  well  as  the  number  of 
differences  taUcen.  The  general  model  for  which  ESTIMATE 
calculates  parameters  is: 

♦ (B)  =*  0^  + 0(B)  a^ 

where  W?  » and  z!  is  the  transformed  value  of  z. . 

ESTIMATE  calculates  estimates  of  the  AR  parameters  for  the 
undifferenced  form  of  the  model,  the  estimated  MA  parameters 
and  MA  constant  term,  the  residual  variance,  autocorrelations 
and  partial  autocorrelations  of  the  residuals,  plots  of 
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these  autocorrelations  and  partial  autocorrelations,  and 
a Chi-square  goodness-of-fit  statistic  for  the  model.  The 
original  series  to  be  modeled  remains  in  FILE  FT02F001 
and  the  residuals  of  the  model  are  stored  in  FILE  FT02F001. 
ESTIMATE  uses  the  IMSL  subroutine  FTMAXL  (a  nonlinear 
gradient  search  algorithm)  to  estimate  the  model  parameters. 

H.  YESTSEAS  PROGRAM 

The  YESTSEAS  progreun  is  used  to  calculate  initial  esti- 
mates of  autoregressive  and  moving  average  parameters  for 
a generalized  seasonal  Box- Jenkins  model.  These  initial 
estimates  are  used  as  inputs  to  the  progreua  WMARQRDT/  which 
then  refines  them  to  determine  the  pareuneter  estimates  that 
minimize  the  sum  of  the  squared  residuals.  The  general 
model  assimed  in  YESTSEAS  is 

<|»(B)  *(B®)7‘^VgZ^  = 9(B)0(B®)a^  . 

The  YESTSEAS  prograim  requires  that  the  seasonal  time  series 
for  which  initial  peurameter  estimates  are  required  is  located 
in  FILE  FT02F001;  it  should  be  tremsformed,  if  desired,  but 
generally  not  differenced,  since  YESTSEAS  performs  the 
required  (seasonal  and/or  non-seasonal)  differencing. 

Terminal  inputs  required  by  the  progreun  include  the  number 
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type  of  pareuneter  to  be  estimated  {that  is,  the  number  of 
each  of  the  seasonal  and  non-seasonal  AR  and  MA  parameters) . 
The  output  of  the  model.  Initial  estimates  for  the  requested 
parameters,  is  printed  out  at  the  terminal.  YESTSEAS  uses 
the  IMSL  subroutines  FTRDIP,  FTMAXL,  FTAUTO  and  FTARPS, 
as  well  as  a modified  version  of  FTMAPS. 

I.  WMARQRDT  PROGRAM 

The  WMARQRDT  program  uses  the  initial  parameter  esti- 
mates calculated  by  YESTSEAS  or  other  means  as  starting  points 
for  calculation  of  parameters  of  a generalized  seasonal 
Box-Jenkins  model  which  minimize  the  sum  of  the  squared 
residuals.  The  general  seasonal  model  assumed  by  WMARQRDT 
can  be  written: 

(()(B)*(B®)v‘^7gZ^  “ ®o  e(B)0(B®)a^. 

Details  of  the  Marquardt-type  non-linear  least  squares 
algorithm  employed  in  this  program  are  provided  in  Chapter 
III  of  this  report.  WMARQRDT  requires  that  the  series  to  be 
modeled  reside  in  FILE  FT02F001;  it  should  be  tremsformed, 
if  desired,  but  generally  not  yet  differenced,  since  the 
program  performs  differencing.  Inputs  at  the  terminal 
include  number  of  seasonal  2md  non-seasonal  differences, 
length  of  the  seasonal  period,  number  and  type  of  parameters 
desired,  and  initial  estimates  of  those  parameters.  The 
program  then  calculates  the  least-squares  parsuneter  estimates. 
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the  steuidard  error  of  these  estimates,  the  MA  constant  term, 
the  residual  variance,  the  autocorrelations  and  partial 
autocorrelations  of  the  residuals,  the  final  model  sum 
of  squared  residuals  and  a Chi>square  goodness-of-fit  sta- 
tistic (the  portmanteau  test) . Additionally,  it  plots 
offline  the  autocorrelations  and  peurtial  autocorrelations  of 
the  residuals.  WMARQRDT  uses  the  IMSL  subroutines  FTRDIF, 
FTAUTO,  LINVIF  and  MDCDFI,  the  SSPLIB  routine  PLOTS,  and 
the  Time  Series  Editor  resident  subroutines  PARSH,  MARQRT, 
SUMSQ,  SWAPS  and  FORMS.  The  calculations  are  performed  in 
double-precision  arithmetic. 

J.  XSUMSQ  PROGRAM 

The  XSUMSQ  program  accpets  an  already  transformed  cuid 
differenced  series  in  FILE  FT02F001  and  general  seasonal 
and  non-seaonal  Box-Jenkins  parameter  values  for  a general 
seasonal  model  from  the  terminal.  The  program  then  calcu- 
lates the  residual  sum  of  squares  for  these  parameters, 
giving  the  user  a feel  for  the  "goodness-of-f it"  the 
parameter  estimates  provide  in  modeling  the  given  time 
series.  XSUMSQ  uses  the  Time  Series  Editor  resident  sub- 
program XSUMSQ.  All  output  is  printed  at  the  terminal. 

The  calculations  eure  performed  in  double-precision  arithmetic. 

K.  FORECAST  PROGRAM 

The  FORECAST  progr2un  uses  the  estimated  Box-Jenkins 
(seasonal  or  non-seasonal)  model  to  compute  forecasts  of 
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the  transformed  and  differenced  time  series.  FORECAST  re- 
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quires  that  the  time  series  to  be  forecasted  reside  in  FILE 
FT02F001  in  transformed  form.  During  execution  the  input 
required  at  the  terminal  includes  the  forecast  origin,  tiie 
numbers  and  estimated  values  of  the  AR  and  MA  parameters 
of  the  model,  the  overall  MA  constant,  the  plot  origin 
index,  the  maximum  forecast  lead  time,  the  order  of  differ- 
encing in  the  model,  and  a level  of  significance  for  the 
forecast  confidence  limits.  The  program  needs  the  treuisfor- 
mation  parameters  from  FILE  FT07F001,  created  when  using 
program  TRANS,  and  performs  the  inverse  transform  to  return 
the  data  to  its  original  form.  FORECAST  also  undiffer- 
ences differenced  data.  The  program  output  includes  the 
forecasts  up  to  the  given  maximum  lead  time,  the  deviations 
from  each  forecast  for  the  confidence  limits,  and  plots  of 
forecasts  and  confidence  limits  plotted  offline.  The 
FORECAST  progreuti  uses  IMSL  subroutine  FTCAST  cuid  the  special  IBM 
SSPLIB  subroutine  UTPLT8.  The  plot  itself  and  output  data 
are  written  onto  FILE  FT08F001  on  the  user's  disc;  since 
this  plot  file  takes  considerable  space,  it  is  usually 
best  to  erase  it  after  execution  of  the  program  is  completed 
and  the  plot  itself  has  been  printed  offline. 

L.  ROOTS  PROGRAM 

The  ROOTS  program  is  used  to  determine  the  roots  of 
the  characteristic  equation  for  a general  ARIMA  model.  The 
roots  are  useful  for  testing  for  stationarity  and  for 
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determining  the  form  of  the  forecast  function.  ROOTS  uses 
the  IMSL  subroutine  ZPOLR  to  calculate  the  roots;  it  requires 
as  input  the  number  of  AR  parameters  in  undifferenced  form,^ 
and  the  actual  values  of  these  parameters,  as  determined 
from  the  ESTIMATE  or  WMARQRDT  programs.  The  roots  are 
printed  out  at  the  terminal.  If  desired,  roots  can  also 
be  obtained  for  the  MA  polynomial,  simply  by  substituting 
the  appropriate  MA  values  for  the  AR  values  as  input. 


M.  GENERATE  AND  SIMULATE  PROGRAMS 

Because  of  their  similarities,  the  GENERATE  and  SIMULATE 
progreuns  are  described  together.  The  GENERATE  progreim  per- 
mits the  user  to  generate  a time  series  from  amy  non-seasonal 
ARIMA  model  he  specifies.  The  user  must  identify  the  model 
and  give  values  for  the  model  parameters  and  starting  con- 
ditions; a random  number  seed  must  also  be  input.  The 
prograun  takes  the  specified  model,  generates  random  noise 
terms,  amd  calculates  as  many  values  of  the  time  series  as 
desired.  The  GENERATE  prograun  can  be  useful  for  purposes 
of  classroom  instruction  in  the  generation  of  a wide  variety 
of  time  series  exaunples  for  model  identification  exercises. 


^Suppose  the  model  was  ARIMA(1,1,0) . The  differenced  form 
would  be  (1-(J),B)V^2  “ + a . Then,  the  undifferenced  form 

would  be  written  (l-(l+<j>,  )B  + (|>,B^)2.  =*  9 + a.  , where  the 

X X u O u ^ 

left-hand  side  was  found  by  multiplying  (l-(|)j^B)  by  (1-B)  » V . 
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It  could  also  be  useful  in  the  diagnostic  phase  of  model 
checkout;  a time  series  could  be  generated  from  the  esti- 
mated model,  and  its  properties  compared  with  those  of 
the  original  series.  If  large  discrepancies  occur  in  this 
comparison;  this  may  be  evidence  that  the  model  is  inade- 
quate. The  model  output,  the  generated  time  series,  is 
written  in  FILE  FT02P001  as  well  as  being  printed  offline. 

The  GENERATE  program  uses  IMSL  subroutine  FTGENl. 

The  SIMULATE  program  provides  the  capability  of  gen- 
erating any  number  of  simulated  time  series.  This  program 
is  useful  for  prediction  of  what  might  happen  in  the  future; 
it  also  demonstrates  that  within  a given  model,  the  time 
series  actually  observed  can  vary  considerably.  This  prograun 
uses  GENERATE,  but  also  requires  as  input  the  actual  series 
in  FILE  FT02F001^  amd  the  number  of  simulated  series  the 
user  wishes  to  generate.  The  series  should  be  already 
transformed  and  differenced.  Additionally,  SIMULATE  allows 
the  user  to  select  v\lues  of  the  original  time  series  as 
starting  values  for  the  simulated  series.  The  output  con- 
sists of  the  simulated  series  printed  at  the  terminal.  The 
SIMULATE  program  uses  IMSL  subroutine  FTGENl. 


N.  HELP  OPTION 

The  HELP  option  is  a modification  to  the  Time  Series 
■ditor  CP/CMS  prograun  itself  that  allows  the  user  to  obtain 


*tf  desired,  the  user  may  choose  to  input  his  own  starting 
mm  and  not  use  the  actual  time  series  values. 
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information  about  Editor  programs  after  the  initial 
introductory  phase  has  been  con^leted. 


0.  SUMMARY 

This  chapter  has  provided  brief  descriptions  of  the 
options  availedsle  for  analysis  of  time  series  data  in  the 
Time  Series  Editor.  The  following  two  chapters  contain 
detailed  examples  of  time  series  analysis  for  both  a non- 
seasonal  and  a seasonal  series;  the  accon^emying  terminal 
sessions  and  computer  output  for  euialysis  of  these  series 
are  contained  in  full  in  Appendices  B and  C.  Appendix  A 
to  this  report  contains  a User's  Guide  for  the  Time  Series 
Editor,  while  Appendix  D contains  a complete  listing  of 
the  Time  Series  Editor  and  resident  programs.  Appendix  E 
contains  listings  of  the  data  sets  used  in  the  analyses  of 
the  next  two  chapters. 
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V.  EXAMPLE  NON-SEASONAL  TIME  SERIES  ANALYSIS 

This  chapter  provides  a description  of  a con^lete  2maly> 
sis  of  a non-seasonal  time  series  using  the  Time  Series 
Editor.  The  time  series  analyzed  is  Series  C (Chemical 
Process  Temperature  Readings),  taken  from  Box  and  Jenkins 
(Ref.  4,  p.  528].  This  time  series  was  selected  since  it 
is  fully  analyzed  and  discussed  in  Reference  4,  so  that  the 
interested  user  can  con5>are  for  himself  the  results  pre- 
sented there  with  the  results  given  by  the  Time  Series 
Editor.  After  reading  this  overview  of  the  analysis,  the 
reader  can  follow  the  actual  user  session  euid  associated 
user  session  Editor  program  output  contained  in  Appendix  B 
to  this  report. 

The  first  step  in  the  analysis  of  series  C is  to  plot 
the  series  using  program  PLOT.  The  plot  reveals  rather  wide 
fluctuations  in  the  series,  but  not  the  sort  of  explosive 
nonstationary  behavior  that  would  render  an  attempt  to  model 
fruitless.  The  plot  also  reveals  that  the  time  series  has 
a large  amount  of  momentum  (movements  of  the  series  tend  to 
resist  changes  of  direction  over  time) . This  is  character- 
istic of  ARIMA  models  with  two  or  more  differences. 

As  a second  step,  the  autocorrelations  (autos)  and 
partial  autocorrelations  (pautos) , mean,  and  Veuriance  of 
the  series  were  estimated  using  AUTO.  The  plot  shows  that 
the  autos  d£unpen  out  slowly  in  a nearly  linear  fashion. 
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This  is  an  indication  that  the  process  is  nonstationary 

and  that  one  or  more  differences  are  required  to  make  it 

stationary.  The  pauto  plot  is  not  informative  when  the 

autos  fail  to  deunpen  out  rapidly. 

As  suggested  by  the  plots  of  the  original  series  2md 

its  auto  and  pautos,  the  program  DIFF  was  used  to  examine 

the  series  with  both  one  and  two  non-seasonal  differences; 

1 2 

i.e.,  the  series  ^ *t  *t  “ ^ *t  derived. 

The  plots  suggest  that  two  differences  might  be  required  to 
achieve  stationarity  in  the  series.  However,  one  might  be 
able  to  get  by  with  a single  difference  and  an  AR  parameter 
near  \inity,  since  the  autos  of  dampen  out  slowly.  Thus, 
two  possible  model  C2uididates  are  suggested: 

''  ' 1 

1.  ARIMA(1,1,0)  : (1  - = ®o  ^t  ' 

eUld 

2.  ARIMA(0,2,2)  : V^z^  * 6^  + (1  - - 02B^)a^  . 

For  purposes  of  estimation,  the  second  model  was  extended 
to  two  moving  average  parcuneters.  Such  overfitting  is  often 
done  to  see  if  the  estimated,  moving  average  parameters  turn 
out  to  be  near  zero,  thus  confirming  the  tentative  identification 
The  next  step  in  the  analysis  process  is  to  calculate 
the  maximum  likelihood  estimates  of  the  model  parameters. 


The  program  ESTIMATE  was  used  to  do  this.  The  estimated 
parameters  for  the  ARIMA(1»1,0)  model  were: 

- 0.8073 

0Q  - -0.006789 
al  - 0.0177633 

A 

The  autos  and  pautos  of  the  residuals  were  calculated 

to  test  the  model.  The  correlations  should  appear  to  be 

estimates  of  a pure  white  noise  process  if  the  model  is 

adequate;  the  (1,1,0)  model  seemed  to  pass  that  test.  The 

2 

chi-square  goodness-of-fit  test  results  in  a x * 28.88  with 
24  df  and  a signficance  level  of  0.2247.  Thus,  there  is 
no  strong  evidence  to  suggest  that  the  (1,1,0)  model  is 
inadequate.  The  model  can  be  written,  in  undifferenced 
form: 


(1  - 1.8073B  + 0.80738“) 2, 


^t  • 


The  9^  is  omitted  here  since  it  is  nearly  zero  in  value. 
The  estimated  parameters  for  the  (0,2,2)  model  were 


0.1378 

0.1296 


e » -0.0026787 
o 

P • 0.0189559 

A 


As  before,  the  correlation  plots  of  the  residuals  fail  to 

suggest  any  inadequacy  of  the  model.  However,  the  chi-square 

2 

goodness-of-fit  test  for  this  model  yielded  ax  * 36.75 
with  23  df  and  a significance  level  of  0.0346.  Therefore, 


i 

i 

I 

I 
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if  the  ARIMA(0,2/2)  model  were  a "correct"  one,  a chi>sguare 
value  as  large  as  36.75  would  occur  by  chance  with  a proba- 
bility of  only  0.0346.  Due  to  its  simplicity  (peursimony 
is  a desirable  feature  in  time  series  models)  and  better 
fit,  the  ARIMA(1,1,0)  model  was  chosen  as  the  "best"  of 
the  two  alternatives  originally  selected. 

The  ARIMA (1,1,0)  model  calculated  was  then  used  to  malce 
forecasts  and  confidence  limits  for  these  forecasts,  using 
program  FORECAST.  The  forecasts  were  made  for  25  periods 
into  the  future,  with  plot  origin  200  and  forecast  origin 
220.  The  length  of  the  series  C itself  is  226;  with  one 
difference,  it  becomes  225.  The  plot  in  Appendix  B shows 
the  forecasted  values,  along  with  the  confidence  limits  at 
level  of  significance  0.10;  the  plot  shows  that  the  width 

' A 

of  the  confidence  limits  increases  substantially  as  the  lag 
gets  large. 

Finally  the  program  ROOTS  was  used  to  calculate  the 
roots  of  the  characteristic  equation  for  the  ARIMA (1,1,0) 

i 

model,  which  can  be  written:  | 

1 - 1.8073B  - 0.8073B^  = 0 j 

•1 

The  roots  calculated  by  ROOTS  were  1.239  and  1.0.  Also,  | 

the  program  SIMULATE  was  used  to  generate  a simulated  series 
from  the  ARIMA (1,1,0)  model,  using  the  last  two  values  of 
the  actual  series  as  starting  values.  The  entire  session 
lasted  about  two  hours,  including  checkout  of  plots  auid 

r- 

')• 

consumed  less  than  a minute  of  CPU  time. 


! I 


87 


t- 

i 


I 


! 


I 


I 


VI.  EXAMPLE  SEASONAL  TIME  SERIES  ANALYSIS 

This  chapter  provides  a description  of  a con^lete  analy- 
sis of  a seasonal  time  series  using  the  Box-JenJcins  tech- 
nique through  the  Time  Series  Editor.  The  time  series 
analyzed  is  Series  G (International  Airline  Passengers: 
Monthly  Totals  (Thousands  of  Passengers)  January  1949  - 
December  1960),  taken  from  Box  euid  Jenkins  [Ref.  4,  p.  531]. 
This  time  series  was  selected  since  it  was  fully  emalyzed 
and  discussed  in  Chapter  9 of  Reference  4;  again,  the 
interested  reader  can  compare  for  himself  the  results  pre- 
sented by  Box  and  Jenkins  with  those  given  by  the  Time 
Series  Editor  programs.  After  reading  this  overview  of  the 
analysis,  the  reader  can  follow  the  actual  user  session  and 
associated  Editor  progr^ua  output  contained  in  Appendix  C 
to  this  report. 

The  plot  shows  both  einnual  seasonality  and  an  increasing 
trend.  This  suggests  that  a seasonal  model  with  differencing 
may  be  needed.  The  autocorrelations  (autos)  and  partial 
autocorrelations  (pautos)  were  then  calculated  and  plotted 
using  AUTO.  The  plot  of  the  autos  shows  very  slow  damping 
out  with  peaks  at  intervals  of  12,  indicating  strong  non- 
stationarity  and  a seasonal  period  of  12  months.  Again, 
here  the  pauto  plot  is  not  informative,  due  to  the  slow 
damping  out  of  the  autos.  As  it  is  sometimes  useful  in 
dealing  with  seasonal  models,  a natural  log  tremsform  was 
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made  on  the  series  G data  using  program  TRANS.  The  plot 
of  the  logged  data  shows  little  chemge  in  non-stationarity, 
and  the  autos  and  pautos  of  the  logged  data  still  indicate 
seasonality  and  non-stationarity. 

Next/  differencing  to  achieve  stationarity  was  tried, 
using  program  DIFF;  both  one  seasonal  and  one  non-seasonal 


difference  of  the  logged  data  was  t2Ucen.  The  autos  and 

,lr,l 


pautos  of  the  data  ^ ^t^  appear  much  inqproved, 

rapidly  danping  out.  They  also  exhibit  sharp  peaks  at  the 
periods  of  seasonality,  indicating  the  need  for  seasonal 
parameterization  in  the  model.  Since  positive  identification 
of  the  model  type  from  the  auto  and  pauto  plots  was  not 
possible,  two  were  postulated  as  close  candidates: 


12 


1.  ARIMA  (0,1,1)  X (0,1,1)  , written: 


V^V^^2(ln  z^)  =60+  (l-BiB)  ' 


emd 


12 


2.  ARIMA  (1,1,1)  X (0,1,1)  , written: 


(1-9,3)  “ ®o  (l-9j^B)  (l-0j^  .,^2®^^^^' 


'1,12' 


The  next  step  in  the  analysis  was  to  calculate  the 
parameter  estimates  for  both  models.  YESTSEAS  was  used  to 
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obteun  initial  parameter  estimates.  For  the  model 


12 


ARINA  (0/1,1)  X (0/1,1)  / YESTSEAS  yielded  starting  values! 


l(s) 


0.390425 


e 


l/12(s) 


0.534397 


For  this  model,  calculation  of  final  parameter  estimates 
in  WMARQROT  using  the  starting  values  above  yielded  the 
following  in  seven  iterations: 


» 0.377152 


®1,12  * 0.572387 


e » 0.000291  ; 

o 


the  chi-square  statistic  for  residual  lack  of  fit  of  the 


model  was  * 29.713571,  with  df  » 38  emd  a probedjility 


of  exceeding  the  x value  of  .829. 


12 


For  the  ARINA  (1,1,1)  X (0,1,1)  model  YESTSEAS  calcu- 
lated starting  values: 


*^1(8)  * 0.112699 


w 


for  this  model,  WMAHQRDT  yielded  the  following  parameter 
values  in  ten  iterations: 

• 0.146667 

» 0.509816 

®1  12  * 0.573388  . 

The  chi-square  statistic  for  residual  lack  of  fit  of  the 

2 . . 
model  was  x “ 28.302543,  with  df  = 37  amd  a probadsxlity 

of  exceeding  the  x value  of  *.847. 

Next,  the  plots  of  autos  and  pautos  of  the  residual  for 

both  models  were  exaunined;  both  models  have  autos  and  pautos 

that  daunpen  out  rapidly,  amd  essentially  appear  to  be 

"white  noise" , with  the  exception  of  a small  peadc  at  auto 

amd  pauto  of  lag  12,  corresponding  to  the  seasonal  period. 

The  chi-square  statistics  indicate  that  both  models  have 

"aUsout"  the  same  goodness  of  fit  to  the  data.  Consequently, 

both  models  were  used  to  make  forecasts  of  the  series. 

Prograun  FORECAST  was  used  to  forecast  series  G into  the 

future  for  30  time  periods,  with  a forecast  origin  of  131; 

the  plot  origin  selected  was  index  100.  As  expected,  both 

models  do  aUsout  the  saune  job  of  forecasting  the  series. 
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The  significance  level  for  the  confidence  limits  for  both 
plots  was  90%. 


5 

i 


As  a final  check  of  both  models,  the  progreun  XSUMSQ 
was  used  to  calculate  the  sum  of  squared  residuals  for 
models.  The  results  were  as  follows: 


1.  For  ARIMA  (0,1,1)  X (0,1,1)^^  , 

S^(9^,  0^^,  ©1^12'  ®a^  “ 0.181929 

2.  For  ARIMA  (1,1,1)  X (0,1,1)^^  , 

9o'®l'®l,12'°a^  * 0.181643  . 


It  can  be  seen  here  that  both  models  are  nearly  equally 
"good",  euid  either  could  be  used  with  a reasonable  amount 
of  confidence  for  forecasting  series  G. 
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There  has  been  an  increasing  need  in  militeury  applies- ■ 
tions,  as  well  as  in  the  general  business  area,  for  a 
computer-aided  time  series  analysis  and  forecasting  capability 
that  is  both  powerful  and  easy  to  use.  Current  military- 
associated  functions  in  which  time  series  techniques  are 
valuable  include  material  inventory  memagement,  recruiting, 
personnel  mcinagement,  budget  analysis,  force  level  projections 
and  perhaps  even  short-term  weather  prediction. 

The  Time  Series  Editor  that  has  been  described  in  this 
report  provides  a unified  collection  of  programs  that  greatly 
facilitate  the  conduct  of  a complete  time  series  analysis 
2md  forecasting  evolution  using  Box- Jenkins  methodology. 

The  routines  guide  the  user  from  the  data  input  stage, 
through  the  data  emalysis  and  parauneter  estimation  stages, 
and  finally  through  the  model  building,  diagnostic  checking 
and  forecasting  stages.  With  its  simple  input  requirements, 
the  Editor  can  be  quickly  mastered  by  even  a beginning 
conqputer  user,  who  has  a basic  understanding  of  Box- Jenkins 
methodology. 

The  Box-Jenkins  methodology  has  been  described  broadly 
in  Chapter  II  and  descriptions  of  the  programs  resident 
in  the  Time  Series  Editor  have  been  given  in  Chapter  IV. 
Exauaple  amalyses  of  both  non-seasonal  and  seasonal  time 
series  are  provided  in  Chapters  V and  VI,  with  complete 


user  sessions  and  output  packages  for  these  emalyses 
included  as  i^pendices  B emd  C.  The  Appendices  also  con- 
tain a Time  Series  Editor  User's  Guide,  as  well  as  con^lete 
program  listings. 

Although  the  Time  Series  Editor  in  its  present  form 
provides  em  excellent  capability  for  interactive  time 
series  analysis  using  the  Box- Jenkins  methodology,  there 
aze  additions  to  the  program  set  that  might  be  made  in  the 
future  which  would  in5)rove  the  Editor's  capability  and 
utility.  It  is  first  recommended  that  the  model  diagnostic 
testing  capability  be  extended  to  include  a periodogram 
analysis,  and/or  other  tests  related  to  the  spectral  analy- 
sis of  time  series  [see  Ref.  4,  p.  294].  Another  useful 
addition  might  be  to  include  in  the  Time  Series  Editor  a 
package  capable  of  iicproving  the  ease  of  entering  datr 
and  simplifying  the  input  error  correction  problem. 
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USER*S  GUIDE  TO  THE  TIME  SERIES  EDITOR 

The  Time  Series  Editor  is  a collection  of  FORTRAN 
programs  driven  by  a control  program  called  TIMESER  EXEC 
in  the  CP/CMS  executive  language  that  has  been  designed 
specifically  for  the  analysis  and  forecasting  of  time 
series  data  using  the  general  Box-Jenkins  methodology. 

This  Guide  contains  information  essential  for  the  user  to 
access  the  Time  Series  Editor,  enter  time  series  data, 
build  a model  emd  obtain  the  desired  output. 

I . DATA  INPUT 

The  Time  Series  Editor  requires  that  time  series  data 
be  entered  into  the  sequential  FORTRAN  input/output  file 
named  FILE  FT02F001.  This  can  be  done  either  online  or 
via  cards  read  offline.  The  following  diagram  illustrates 
the  proper  card  deck  arrangement.  When  your  deck  is  ready 
for  input,  give  it  to  a computer  center  system  operator 
to  direct  it  to  your  disc  space. 

If  the  user  has  a data  deck  already  punched  up  in  a 
format  other  than  5F15.6,  he  may  enter  the  series  as  above 
into  FILE  FT03F001  (without  the  length  of  series  card)  emd 
use  the  Editor  program  ZFORMAT  to  transform  it  into  the 
FILE  FT02F001  in  the  proper  format,  without  destroying  his 
original  file.  This  program  is  described  in  the  next  section. 
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II.  TABLE  OF  OPTIONS 

This  Table  provides  the  user  with  the  basic  inforxnation 
necessairy  to  understand  the  data  requirements,  functions 
and  output  options  for  each  program  in  the  Editor. 


Files;  Files:  (1)  al 


Files 
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III.  THE  BASIC  USER  SESSION 

To  use  the  Time  Series  Editor,  the  user  must  log  into 
CMS,  get  into  CP,  link  to  the  disc  storage  area  where  the 
Time  Series  Editor  resides,  reimplement  CMS,  log  into  the 
general  user  and  Time  Series  Editor  disc  areas,  and  enter 
the  TIMESER  routine.  This  section  will  provide  explicit 
guidelines  to  enedjle  the  user  to  perform  the  above  steps 
on  the  NFS  CP/CMS  system.  Commands  marked  with  an  asterisk 
(*)  aure  those  actually  entered  on  the  terminal  by  the  user 
(the  asterisk  itself  is  omitted) . Those  without  an  aster- 
isk euid  those  written  in  all  capital  letters  are  system 
responses  at  the  terminal.  Numbered  sentences  are  comments, 
which  will  not  appear  during  an  actual  user  session.  The 
instructions  and  system  responses  assure  the  user  is  on  am 
IBM  2741  Input/Output  Terminal.  Some  minor  modifications 
may  be  necessary  if  other  terminals  are  used. 


Turn  the  terminal  on,  depress  the  BREAK  key,  and  wait 
for  the  system  to  respond: 

CP-67  online  xd.65  gsyosu 

Depress  the  ATTN  key.  The  roll  bar  will  advance  and 
the  keyboard  will  unlock.  Then  enter: 

* login  aaaapbb  450k 

aaaa  is  the  user's  identification  number,  amd  nn  is  the 
terminal  number  (usually  written  on  the  terminal) . 

For  example,  if  the  user's  ID  number  is  1621  amd  the 
terminal  numbexf  is  44,  the  input  would  be: 
login  1621p44  450k.  The  addition  of  450  k to  the  normal 
login  commamd  is  necessairy  to  execute  the  program 
WMARQRDT  in  the  Editor;  for  users  not  plamning  to 
execute  this  program  during  a session,  this  addition 
is  not  necessary. 
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The  system  will  respond  with  the  statement: 

ENTER  PASSWORD: 

The  user  then  enters  his  password,  or  the  general  users 
password  npg; 

*pas sword 

The  system  will  then  respond: 

ENTER  4-DIGIT  PROJECT  NUMBER  FOLLOWED  BY  4 -CHARACTER 
COST  CENTER  CODE: 

The  user  then  enters: 

♦gggghhhh 

gggg  is  the  assigned  project  number,  and  hhhh  is  the 
user's  section  designator  or  the  faculty  code. 

The  system  will  respond  with  the  message  of  the  day, 
such  as : 

CP/CMS  HOURS  ...  0930=2200 (MON-THURS)  ...  0930-1800 (FRI) 
•OUTPUT  RETAINED  5 DAYS 
Cms  Version  3.25 

At  this  point,  the  user  is  in  CMS.  He  must  then  get  into 
CP;  this  can  be  done  by  hitting  the  ATTN  key.  The  system 
will  then  respond: 


The  user  must  then  link  to  the  TIME  SERIES  EDITOR; 
this  is  accomplished  by  entering: 

♦link  2069p  191  192 

The  system  will  respond  with: 

ENTER  PASSWORD: 

The  password  (read  only)  to  enter  the  Editor  is: 
♦timser 

The  system  then  responds: 

SET  TO  READ  ONLY 


15.  The  user  now  implements  CMS  by: 

*ipl  cms 

16.  The  system  will  respond: 

CMS  Version  3.25 

17.  Now  the  user  must  log  into  both  the  general  user  and 
the  Time  Series  Editor  area,  by  entering: 

♦login  191 

18.  The  system  will  respond  with  a message  such  as: 


19.  The  user  then  enters  the  command: 

♦login  192  t,p 

20.  The  system  will  respond: 

T (192)  R/0 

R; 

21.  The  user  ceui  then  enter  the  Time  Series  Editor 
(guided  version)  by  entering  the  corarMuid; 

♦timeser 

22.  The  system  will  respond: 

EACH  2 SECONDS  EXECUTION  TIME  IS  INDICATED  BY  ♦ 

YOU  HAVE  ENTERED  THE  TIME  SERIES  EDITOR 

PLEASE  RESPOND  TO  EACH  QUERY  WITH  AN  INPUT  AT  THE  TERMINAL. 
ENTER  ONLY  THE  FIRST  LETTER  FOR  A WORD  RESPONSE. 

ENTER  NUMERICAL  VALUES  VIA  FORTRAN  FORMAT.. 

TYPE  INTEGER  VALUES  (RIGHT  JUSTIFIED)  FOR  NAMES  STARTING 
WITH  I THROUGH  N.  TYPE  FLOATING  VALUES  WITH  DECIMAL  FOR 
ALL  OTHERS. 

DO  YOU  WANT  A LIST  OF  THE  OPTIONS? 

23.  The  user  is  then  on  his  own,  guided  by  the  Exec  routine. 
See  the  notes  that  appear  at  the  end  of  this  guide  for 
additional  information.  Eventually  the  user  will  be 
asked: 

DO  YOU  WANT  TO  TRY  AGAIN? 


24.  If  a yes  response  Is  given,  zmother  sequence  will  begin; 
if  the  response  is  no,  the  user  will  be  taken  out  of 
the  Time  Series  Editor  environment  and  returned  to  CMS. 
The  system  response  will  be: 

CONTROL  RETURNED  TO  CMS 
R; 

25.  The  user  can  then  log  out  of  CMS  by  typing: 

*cp  logout 

26.  The  system  will  respond  with: 

CONNECT-  00:08:02  VIRTCPU-  000:07.98  TOTCPU-  000.10.94 
LOGOUT  AT  14.22.04  on  10/16/78 

27.  The  user  should  then  turn  off  his  terminal  and  tear  off 
the  output  from  his  session. 

The  more  experienced  user  can  dispense  with  the 
"welcome  aboard"  section  of  the  Time  Series  Editor  and  get 
right  down  to  business  by  using  the  shortened  version  of  the 
Editor.  This  shortened  version  may  be  entered  by  linking  in 
the  normal  way,  and  then  entering  the  Editor  by  typing  the 
COMMAND 


*timeser  s (asterisk  omitted) . 


The  system  will  immediately  respond: 


ENTER  LETTER  FOR  OPTION  YOU  WANT. 


The  session  inside  the  Editor  then  begins. 


A brief  sample  user  session  is  given  below;  it  includes 
copies  of  the  offline  output  generated  during  the  session. 


repeat  login  nucfpnff 

login  1621p44  450k 
ENTER  PASSWORD < 

ENTER  4-DIGIT,  PROJECT  NUMBER  FOLLOWED  BY  4-CHARACTBR  COST  CENTER  COOEi 
0444rl72 

READY  AT  17.24.38  ON  08/16/78 
CMS  Version  3.25 


P (191)1  29  FILBSi  241  REC  IN  OSEr  55  LEFT  (oC  296),  81%  FULL  (2  CYL) 
R» 

cp  q C 

FILESs-  NO  RDR,  NO  PRT,  NO  PUN 

Rl 

CP 

link  2069P  191  192 
ENTER  PASSWORD t 

SET  TO  READ  ONLY 

ipl  erne 

CMS  Version  3.25 


login  191  t 

R»  ! 

login  192  t,p 

T (192)  R/0  1 

i 


timeser 

EACH  2 SECONDS  EXECUTION  TINE 


IS  INDICATED 


BY 


YOU  HAVE  ENTERED  THE  TIME  SERIES  EDITOR. 


PLEASE  RESPOND  TO  EACH  QUERY  WITH  AN  INPUT  AT  THE  TERMINAL'. 
ENTER  ONLY  THE  FIRST  LETTER  FOR  A WORD  RESPONSE. 

ENTER  NUMERICAL  VALUES  VIA  FORTRAN  FORMAT. 


TYPE  INTEGER 
WITN  I THRU  N 


VALUES  (RIGHT  JUSTIFIED)  FOR  NAMES  STARTING  ] 

. TYPE  FLOATING  VALUES  WITH  DECIMAL  FOR  ALL  OTHERS. 


DO  YOU  WANT  A LIST  OF  THE  OPTIONS? 

y 


OPTION  DESCRIPTION 

GENERATE  -GENERATE  ANY  ARINA  TIME  SERIES 

AUTO  - — “ — —CALCULATE  AUTOCORRELATIONS,  PAUTOS,  MEAN  AND  VARIANCE 
PLOT  ———PLOT  A TIME  SERIES 

ESTIMATE  —CALCULATE  MAX  LIKELIHOOD  ESTIMATES  OP  ARMA  PARAMETERS 
DIfp DIFFERENCE  A TIME  SERIES 

FORECAST  FORECAST  FUTURE  VALUES,  CONSTRUCT  CONFIDENCE  INTERVALS 

TRANS •TRANSFORMS  VALUES  OF  A TIME  SERIES 

ROOTS  ——DETERMINES  ROOTS  OF  ARINA  CHARACTERISTIC  EQUATION 
ZFORMAT  ALTER  DATA  FILE  TO  FORMAT  5F15.6 

CNSWORK  ———PERFORM  CP/CNS  COMMANDS  IN  TIMBSEi*  EXEC 

SIMULATE  — SIMULATE  NONSBASONAL  TINE  SERIES 

YBSTSEAS  ——CALCULATE  INITIAL  SEASONAL  PARAMr’BRS 

WMAROROT  -MARQUARDT  SOLUTION  FOR  PARAMETER  ESTIMATES 

XSUMSQ  ——CALCULATE  SOM  OF  SQUARES  FOR  ARBITRARY  PARAMETERS 


WOULD  YOU  LIKE  MORE  INFO? 

y 


I 


110 


ENTER  OPTION  TOO  WANT  INFO  ABOUT 

a 


AUTO  THIS  PROGRAM  CALCULATES  AUTOCORRELATIONS,  PARTIAL 

AUTOCORRELATIONS,  THE  MEAN  AND  TUB  VARIANCE  FOR  A GIVEN  TINE  SERIES 
nHICH  MOST  RESIDE  IN  FILE  FT02F001.  THE  PROGRAM  OSES 
FTAUTO  IN  THE  IMSL  LIBRARY.  THE  AUTOCORRELATIONS  AND  PAUTOS  CAM  BS 
PLOTTED  OFFLINE. 


DO  YOU  WANT  INFO  ABOUT  ANOTHER  OPTION? 

n 

DO  YOU  WANT  TO  TRY  A SESSION? 

y 

ENTER  LETTER  FOR  OPTION  YOU  WANT, 

c 

ENTER  DESIRED  CP/CMS  COMMANDS,  ONE  PER  LINE. 
WHEN  FINISHED  TYPEi  tGOTO  -QUES 
liste  * ft02£00X 


FILENAME 

FILETYPE 

MODE 

NO. REC. 

DATE 

SBRG 

FT02F001 

PI 

3 

9/16 

SERC 

FT02F001 

PI 

5 

9/16 

LNSERG 

FT02F001 

PI 

3 

9/16 

FILE 

FT02F001 

PI 

1 

9/16 

cease  file  ft02I001 

alter  sere  £t02£001  pi  file  £t02£001  pi 
stat 

P (191):  28  FILES}  240  REC  IN  OSS,  56  LEFT  (o£  296),  81%  FULL  (2  CYL) 
iigoto  -ques 

DO  YOU  WANT  TO  GO  AGAIN? 

y 


ENTER  LETTER  FOR  OPTION  YOU  WANT. 

a 

IS  YOUR  DATA  IN  FILE  PT02F001? 

y 

EXECUTION  BEGINS... 

AUTOCORRELATIONS 

0.978  0.944  0.902  0.854  0.802  0.748  0.692  0.635  0.579  0.923 

0.468  0.413  0.359  0.305  0.253  0.201  0^150  0.098  0.047  -0.003 

-0.052  -0.101  -0.151  -0.200  -0.248 

PARTIAL  AtrroCORRELATIONS 

0.978  -0.260  -0.157  -0.093  -0.058  -0.045  -0.012  -0.038  -0.022  -0.010 

-0.036  -0.041  -0.038  -0.024  -0.037  -0.027  -0.032  -0.070  -0.048  -0.024 

-0.034  -0.061  -0.079  -0.048  -0.037 


MEAN-  22.9739 


VARIANCE  - 4.22273 


ENTER  TITLE  FOR  PLOTS, 
autos  and  pautos  for  series  c data 
YOUR  AUTO  AND  PAOTO  PLOTS  HAVE  BEEN  PRINTED  OFFLINE. 
PICE  OP  IN  ROOM  1140  UNDER  YOUR  USER  ID  NUMBER. 

DO  YOU  WANT  TO  GO  AGAIM? 

y 

ENTER  LETTER  FOR  OPTION  YOU  WANT. 

C 

ENTER  DESIRED  CP/CMS  COMMANDS,  ONE  PER  LINS. 

WHEN  FINISHED  TYPE:  tCOTO  -QUES 
Offline  print  file  ft02£001 
alter  file  ft02f001  pi  sere  ft02f001  pi 
agoto  -ques 

DO  YOU  WANT  TO  GO  AGAIN? 
n 

CONTROL  RETURNED  TO  CMS 
R} 
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V.  PROBLEM  CONTROL  NOTES 

This  section  will  cover  corrective  measures  that  cam  be 
taken  when  things  fail  to  go  as  expected  while  in  the 
TIMESER  environment. 

a.  A typing  error  in  the  CMS  environment  can  be 
corrected  by  typing  the  @ character  as  mamy  times  as  is 
required  to  back  up  and  then  type  the  correct  values.  For 
examqpler  if  the  user  typed  times re/  the  user  could  correct 
the  mistake  by  typing  two  @ signs,  followed  by  the  correct 
spelling  er,  as  follows:  timesre88re.  An  entire  line  can 

be  deleted  by  typing  the  ^ character  (or  [ on  some  terminals) 

b.  When  working  with  TIMESER  executive  prograuns,  the 
user  should  exercise  care  before  hitting  the  return  key  (or 
control s on  some  terminals) . If  an  input  value  is  required 
and  the  return  key  is  hit  before  the  proper  response  is 
entered,  the  user  will  be  likely  to  get  thrown  out  of  the 
editor  amd  have  to  begin  again.  In  most  cases,  errors  can 
be  corrected  only  before  the  return  key  is  struck  (in  some 
programs  you  get  a second  chamce  for  input) . 

c.  Particular  caure  should  be  taUcen  for  integer  value 
input,  which  must  be  right  justified  in  the  format  field. 

The  editor  will  advise  the  user  in  all  cases  where  the 
Integer  format  is  other  tham  II. 

d.  If  for  amy  reason  the  user  finds  himself  in  a debug 
or  error  condition  (caused  by  erroneous  data,  or  a "blowup" 
in  one  of  the  non-linear  optimization  routines  usually 


116 


caused  by  very  poor  initial  input  values) , the  following 
procedure  will  get  the  user  back  into  the  normal  TIMESER 
environment: 

(1)  depress  the  ATTN  key  twice;  this  gets 
the  user  into  CP;  hit  the  ATTN  key  again, 
emd  then  type  kx  to  kill  the  execution. 

(2)  re-ipl  CMS,  and  login  191  and  then  login 
192  t,p. 

(3)  then  the  user  can  type  TIMESER  or  TIMES^  S, 
amd  return  to  the  TIMESER  environment 

On  the  next  page  is  a sample  user  session  where  an  error 

causing  a debug  condition  has  occurred. 

Here  the  user  executed  a prograun  that  required  data 
in  FILE  FT02F001,  and  the  FILE  did  not  exist.  As  the 
example  shows,  recovery  is  quick.  Simply  hit  the  break 
button  to  get  into  CP,  login  191,  login  192  t,p,  and  then 
type  timeser  s to  return  immediately  to  the  Editor  environment. 


ENTER  DESIRED  CP/CMS  COMMANDS,  ONE  PER  LINE 
WHEN  FINISHED  TYPE;  &GOTO  -QUES 
alter  sere  ft02f001  pi  file  ft02f001  pi 
&goto  -ques 

DO  YOU  WANT  TO  GO  AGAIN? 


APPENDIX  B 


KEYBOARD  SESSION  AND  OUTPUT  FROM  ANALYSIS  OF 
TIME  SERIES  C (NON-SEASONAL) 

This  appendix  contains  the  complete  terminal  keyboard 
session  and  all  offline  output  from  the  euialysis  of  Box 
and  Jenkins'  [Ref.  4]  time  series  C (Chemical  process 
temperature  readings)  using  the  Time  Series  Editor.  This 
appendix  is  intended  to  supplement  the  discussion  of  the 
analysis  of  series  C discussed  in  Chapter  V of  this  report. 

The  terminal  listings  are  presented  first,  followed  by 
the  conqputer  offline  output. 

kE  i 


repeat  lo9in  nurfpnf 

login  l(21p44  400k 
ENTER  PASSMOROt 

ENTER  4>DXGIT  PROJECT  NONBBR  POLLONBD  BY  4-CaARACTBR  COST  CENTER  COOBi 
0444rl72 

READY  AT  14.23.02  ON  09/16/78 
CMS  vereieif  3.25 


1 listf  * 

£t02f001 

1 FILENAME 

FILETYPE 

NODE 

NO.RBC. 

DATE 

r FILE 

FT02F001 

PI 

3 

9/15 

SERG 

PT02F001 

PI 

3 

9/16 

SERC 

PT02f001 

PI 

5 

9/16 

eceee  file  £t02£001 
R> 

alter  sere  £t02£001  pi  file  £t02£001  pi 
Ri 

CP 

link  2069p  191  192 
ENTER  PASSWORD: 

SET  TO  READ  ONLY 

ipl  cm 

CHS  Version  3.25 

login  191 
Ri 

login  192  t.p 
T (192)  R/0 

R;  

^er 

..aCH  2 SECONDS  EXECUTION  TINE  IS  INDICATED  BY  *. 

YOU  HAVE  ENTERED  THE  TIME  SERIES  EDITOR. 

PLEASE  RESPOND  TO  EACH  QUERY  WITH  AN  INPUT  AT  THE  TERHINAL. 

ENTER  ONLY  THE  FIRST  LETTER  FOR  A WORD  RESPONSE. 

ENTER  NUMERICAL  VALUES  VIA  FORTRAN  FORMAT. 

TYPE  INTEGER  VALUES  (RIGHT  JUSTIFIED)  FOR  NAMES  STARTING 

WITH  I THRU  N.  TYPE  FLOATING  VALUES  WITH  DECIMAL  FOR  ALL  OTHERS. 

DO  YOU  WANT  A LIST  OF  THE  OPTIONS? 

y 

OPTION  DESCRIPTION 

GENERATE  -GENERATE  ANY  ARIMA  TIME  SERIES 

AUTO CALCULATE  AUTOCORRELATIONS,  PAUTOS,  MEAN  AND  VARIANCE 

PLOT PLOT  A TIME  SERIES 

ESTIMATE  -CALCULATE  MAX  LIKELIHOOD  ESTIMATES  OF  ARHA  PARAMETERS 

DIPF  — — DIFFERENCE  A TIME  SERIES 

FORECAST  — — — FORECAST  FUTURE  VALUES,  CONSTRUCT  CONFIDENCE  INTERVALS 

TRANS  TRANSFORMS  VALUES  OF  A TIME  SERIES 

ROOTS  — — — DETERMINES  ROOTS  OF  ARINA  CHARACTERISTIC  EQUATION 

ZFORMAT  ALTER  DATA  FILE  TO  FORMAT  5F15.6 

CMSWORR  PERFORM  CP/CHS  COMMANDS  IN  TIMESER  EXEC 

SIMULATE  — SIMULATE  NONSEASONAL  TIME  SERIES 

YESTSBAS  — >-CALCULATB  INITIAL  SEASONAL  PARAMETERS 

HMARQROT  MARQUARDT  SOLUTION  FOR  PARAMETER  ESTIMATES 

XSUMSQ  — CALCULATE  SON  OP  SQUARES  FOR  ARBITRARY  PARAMETERS 


WOULD  YOU  LIRE  MORE  INFO? 

y 


'H 


II I Ml 
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ENTER  OPTION  TOO  NANT  INPO  ABOUT 

a 


AUTO  THIS  PROGRAM  CAI.COLATBS  AUTOCORRELATIONS,  PARTIAL 

AUTOCORRELATIONS,  THE  MEAN  AND  THE  VARIANCE  POR  A GIVEN  TIME  SERIES 
WHICH  MUST  RESIDE  IN  PILE  PT02P001.  THE  PROGRAM  OSES 
PTAUTO  IN  THE  INSL  LIBRARY.  THE  AUTOCORRELATIONS  AND  PAOTOS  CAN  BE 
PLOTTED  OPPLINE. 


DO  YOU  NANT  INPO  ABOUT  ANOTHER  OPTION? 

y • 

ENTER  OPTION  YOU  WANT  INPO  ABOUT, 
d 

DIPP  THIS  PROGRAM  TARES  A SERIES  IN  PILE  PT02P001,  ASKS  IP 

IT  IS  SEASONAL,  AND  THEN  PERFORMS  NONSEASONAL  AND/OR  SEASONAL 
OIPPBRENCING,  AS  REQUESTED.  THE  ORIGINAL  DATA  REMAINS  IN  PILE 
PT02P001,  WHILE  THE  DIPPERENCBD  SERIES  IS  POT  IN  PILE  PT03P001. 


DO  YOU  WANT  INPO  ABOUT  ANOTHER  OPTION? 

n 

00  YOU  WANT  TO  TRY  A SESSION? 

Y 

ENTER  LETTER  POR  OPTION  YOU  WANT. 

a 

IS  YOUR  DATA  IN  PILE  PT02P001? 

y 

EXECUTION  BEGINS... 

AUTOCORRELAT IONS 


0.978  0.944 

0.468  0.413 

-0.052  -0.101 


0.902  0.854 
0.359  0.305 
0.151  -0.200 


0.802  0.748 

0.253  0.201 

-0.248 


0.892  0.635 

0.150  0.098 


0.579  0.523 

0.047  -0.003 


PARTIAL  AUTOCORRELATIONS 

0.978  -0.260  -0.157  -0.093  -0.058  -0.045  -0.012  -0.038  -0.022  -0.010 
-0.036  -0.041  -0.038  -0.024  -0.037  -0.027  -0.032  -0.070  -0.048  -0.024 
-0.034  -0.061  -0.079  -0.048  -0.037 


MEAN*  22.9739  VARIANCE  - 4.22273 

ENTER  TITLE  POR  PLOTS. 

autos  and  pautos  of  secies  c data  / undiffecenced 
YOUR  AUTO  AND  PAUTO  PLOTS  HAVE  BEEN  PRINTED  OFFLINE. 
PICK  OP  IN  ROOM  1140  UNDER  YOUR  USER  ID  NUMBER. 

*DO  YOO  WANT  TO  GO  AGAIN? 

y 


BNT£R  LETTER  FOR  OPTION  YOU  WANT, 
d • 

IS  YOOR  DATA  IN  PILE  PT02P0017 

y ^ 

EXECOTION  RB61NS... 

IS  YOUR  TINE  SERIES  SEASONAL? 

n 

ENTER  NOMSBR  OP  NOMSBASONAL  DIFFERENCES. 

1 

•DO  YOU  WANT  TO  PLOT  AUTO  AND  PAUTO  OP  TRANSPORHBD  DATA? 

y 

EXECUTION  BEGINS... 

AUTOCORRELATIONS 


0.805 

0.653  0.526 

0.442 

0.380 

0.318 

0.262 

0.186 

0.139 

0.144 

0.097 

0.094  0.074 

0.073 

0.070 

0.072 

0.089 

0.048 

0.041 

0.040 

0.044 

0.048  0.000 

-0.052 

-0.097 

PARTIAL  autocorrelations 

0.805 

0.010  -0.007 

0.051 

0.027 

-0.019 

-0.013 

-0.080 

0.020 

0.117 

-0.137 

0.094  -0.027 

0.034 

0.006 

0.012 

0.043 

-0.121 

0.060 

-0.005 

0.038 

-0.022  -0.116 

-0.076 

-0.008 

1 


MEAN— .346667E-01  VARIANCE  « 0.S31982E-01 

ENTER  TITLE  POR  PLOTS. 

autos  and  pautos  of  secies  c data  / one  nonseasonal  difference 

YOOR  AUTO  AND  PAUTO  PLOTS  HAVE  BEEN  PRINTED  OFFLINE. 

PICK  OP  IN  ROOM  1140  UNDER  YOUR  USER  ID  NUMBER. 

DO  YOU  WANT  TO  GO  AGAIN? 

y 

ENTER  LETTER  POR  OPTION  YOU  WANT, 
d 

IS  YOUR  DATA  IN  PILE  PT02F001? 

y 

EXECUTION  BEGINS... 

IS  YOUR  TINE  SERIES  SEASONAL? 

n 

ENTER  NUMBER  OP  NONSEASONAL  DIFFERENCES. 

2 

DO  YOU  WANT  TO  PLOT  AUTO  AND  PAUTO  OP  TRANSFORMED  DATA? 

y 

EXECUTION  BEGINS... 

AUTOCORRELATIONS 


-0.079 

-0.122 

-0.007 


-0.06S 

0.072 

0.135 


-0.122 

-0.077 

0.014 


-0.063 

0.029 

-0.013 


0.013 

-0.011 

-0.135 


-0.018 

-0.058 


0.049  -0.052  -0.124  0.122 

0.171  -0.101  -0.013  -0.020 


PARTIAL  AUTOCORRELATIONS 

-0.079  -0.072  -0.135  -0.094  -0.022  -0.051  0.022  -0.060  -0.145  0.094 
-0.143  0.023  -0.089  -0.005  -0.039-0.074  0.130  -0.105  -0.014  -0.053 
0.004  0.081  0.056  -0.041  -0.099 


MEAN— .267866E-02  VARIANCE  - 0.198143B-01 


ENTER  TITLE  POR  PLOTS. 

autos  and  pautos  of  series  c data  / two  nonseasonal  differences 

YOUR  AUTO  AND  PAUTO  PLOTS  HAVE  BEEN  PRINTED  OFFLINE. 

PICK  OP  IN  ROOM  1140  UNDER  YOUR  USER  ID  NUMBER. 


i 

1 

I 
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1 


i i 


DO  you  WANT  TO  GO  AGAIN? 

y 

BNTBR  bBTfBit  rOR  OPTION  TOO  WANT. 

C 

BNTBR  DESIRED  CP/CMS  COMMANDS.  ONB  PER  LIMB. 

WHEN  PIN^SHBO  TTPBt^  (GOTO  -QUES 
«lt«c  Cil*  ftOlfOQI  pX  Mce  ftozeool  pi 
«lt«t  file  ftOSTOOl  pi  file  ft02f001  pi 
AGPto  -qaas 

DO  YOU  WANT  TO  GO  AGAIN? 

y 

ENTER  LETTER  FOR  OPTION  YOU  WANT. 

IS  YOUR  DATA  IN  PILE  PT02F001? 

y 

BX8CUTZON  BEGINS* •• 

ENTER  NOMBBR  OPM  PARAMETERS  (EXCLUDE  DIFFERENCES). 

0 

BNTBR  NUMBER  OF  MA  PARAMETERS. 

2 

LENGTH  OF  TIME  SERIES  - 224 
0 AR  PARAMETERS  2 MA  PARAMETERS 

MA  PARAMETERS  ARE: 

THBTA(l)-  0.1378 


THBTA{2)-  0.1296 
MA  CONSTANT-  -0.26787E-02 


WHITE  NOISE  VARIANCE- 
36.75  OF 


0.189S59B*01 


CHI  SQUARE  RESIDUAL  LACK  OF  FIT  VALUE  - 36.75  OF  - 23 

SIGNIFICANCE  - 0.0346 

DO  YOU  WANT  TO  PLOT  AUTO  AND  PAUTO  OF  RESIDUALS? 

y 

^EXECUTION  BEGINS... 

AUTOCORRELATIONS 

0.019  0.034  -0.134  -0.093  -0.012  -0.043  0.005  -0.063  -0.125  0.088 
-0.142  0.069  -0.098  0.025  0.004  -0.055  0.157  -0.090  0.005  -0.006 
0.009  0.133  0.005  -0.015  -0.133 


PARTIAL  AUTOCORRELATIONS 

0.019  0.033  -0.135  -0.090  0.001  -0.055 

-0.166  0.017  -0.104  -0.013  -0.025  -0.089 
0.024  0.069  0.024  -0.069  -0.102 


-0.019  -0.071  -0.143  0.087 

0.123  -0.123  -0.021  -0.025 


MBAM-0.236684E-02 


VARIANCE  - 0.189503B-01 


ENTER  title  FOR  PLOTS. 

autos  Slid  pautos  oC  casiduals  / ariu(0,2,2)  Model  ot  sexles  e data 
YOUR  AUTO  AND  PAUTO  PLOTS  HAVE  BEEN  PRINTED  OPPLXNB. 

PICK  OP  IN  ROOM  1140  UNDER  YOUR  USER  ZD  NUMBER. 

DO  YOU  WANT  TO  GO  AGAIN? 

y 


123 


I 


ENTER  LETTER  FOR  OPTION  TOO  WANT, 
c 

ENTER  0E8XRS0  CP/CN8  COHMANOS,  ONI  PCS  LIU. 

WHEN  PZNXSHSD  TYPES  MSOTO  -QOE* 

•CAM  fil*  £t02C001 

P^(191)i  27  PZLBEl  2CS  REC  ZE  OSE,  21  LEFT  (of  29(),  90%  POLL  (2  CYL) 
•CAM  CtOECOOl 

Aie*c  ••re  £t02C001  pi  Cil*  £t02£001  pi 
(goto  -qUAS 

00  YOO  WANT  TO  60  A6AZN? 

y 

ENTER  LETTER  POR  OPTION  YOU  NANT. 
d 

•ZS  YOUR  DATA  ZN  PZLE  PT02P0017 

y 

I BXBCUTZON  BEGINS... 

IS  YOUR  TZHE  SERIES  SEASONAL? 
n 

I ENTER  NUMBER  OP  NONSEASONAL  OIPPERENCES. 

’ 1 

I 00  YOU  WANT  TO  PLOT  AUTO  AMO  PAOTO  OP  TRAMSPORNEO  DATA? 

* n 

00  YOO  NANT  TO  GO  AGAIN? 

y 

ENTER  LETTER  POR  OPTION  YOU  WANT. 

: . c 

•ENTER  OBSIRBO  CP/CMS  COHNANOS.  ONE  PER  LINE. 

WBBN  PINISHEO  TYPE:  AGOTO  -QOES 
Altec  file  £t02f001  pi  SAcc  £t02£001  pi 
AltAC  CilA  ft03C001  pi  file  £t02£001  pi 
I Aqoto  -^ucs 

DO  YOU  NANT  TO  GO  AGAIN? 


ENTER  LETTER  POR  OPTION  YOO  WANT. 

A 

IS  YOUR  DATA  IN  PILE  PT02P001? 

I y 

EXECUTION  BEGINS... 

i ENTER  NUMBER  OP  AR  PARAMETERS  (EXCLUDE  OIPPERENCES). 

- 1 

P ENTER  NUMBER  OP  NA  PARAMETERS. 

0 

• LENGTH  OP  TINE  SERIES  - 225 
1 AR  PARAMETERS  0 MA  PARAMETERS 

AR  PARAMETERS  AREs 
PHZ(1)>  0.9073 

NA  CONSTANT-  -0.697990-02  WHITE  NOISE  VARIANCE-  0.177633B-01 

CBZ  SQUARE  RESIDUAL  LACE  OP  PIT  VALUE  - 29.99  OP  - 24 
SI6MIPICANCE  - 9.2247  ’ ■ 

DO  YOU  WANT  TO  PLOT  AUTO  AND  PAOTO  OP  RESIDUALS? 

I EXECUTION  BEGINS... 

AUTOCORRELATIONS 

0.019  0.014  -0.049  -0.005  0.091  0.024  0.079  -0.019  -0.097  0.130 
-0.093  0.094  -0.055  0.043  0.007  -0.040  0.199  -0.079  -0.001  -0.010 
0.000  0.129  0.014  -0.013  -0.127 


PARTIAL  AUTOCORRELATIONS 


HBAM— .2S2S17B-02  VARIANCE  - 0.177569B-01 

ENTER  TITLB  FOR  PLOTS.  . i ^ 

autos  and  pautoa  of  rasiduala  / afi«a(l»l»0)  nodal  of  aaciaa  e data 

YOUR  AUTO  ANO  PAOTO  PLOTS  HAVB  BBBN  PRINTBO  OPPLXNB. 

PICB  UP  IN  ROOM  X140  UNDER  YOUR  USER  ID  NUHBBR. 

DO  YOU  NANT  TO  GO  AGAIN? 


ENTER  LBTTBR  FOR  OPTION  YOU  WANT 


ENTER  DESIRED  CP/CMS  COMHANDS,  ONE  PBR  LINE 
WHEN  FINISHED  TYPSi  (GOTO  -QUBS 
eras*  file  ft02f001 

alter  aerc  ft02f001  pi  file  ft02f001  pi 
(goto  *^ues 

DO  YOU  WANT  TO  GO  AGAIN? 


ENTER  LETTER  FOR  OPTION  YOU  WANT 


IS  YOUR  DATA  IN  FILE  FT02F0017 


EXECUTION  BEGINS 


IS  YOUR  SERIES  SEASONAL? 


ENTER  HONSEASONAL  DATA  AS  REQUESTED t 


ENTER  NUMBER  OF  NON-SBASONAL  DIFFERENCES 


ENTER  HUMBER  OF  MON-SEASONAL  AR  PARAMETERS 


ENTER  NUMBER  OF  NON-SBASONAL  MA  PARAMETERS 


NOW  INPUT  YOUR  PARAMETER  ESTIMATBSt 


ENTER  MA  CONSTANT  TERM,  THETAO. 

.00667 

ENTER  NON-SBASONAL  AR  PARAMETER  PHI(l) 


ENTER  NAXXNUN  FORECAST  LEAD  TINE  VIA  FORMAT  12. 

25 

ENTER  INDEX  FOR  PLOT  ORIGIN  VIA  FORMAT  13. 

200 

ENTER  INDEX  FOR  FORECAST  ORIGIN  VIA  FORMAT  13. 

220 

ENTER  SIGNIFICANCE  LEVEL  FOR  CONFIDENCE  INTERVALS 

0.10 

WAS  YOqR  DATA  TRANSFORMED  IN  THE  TRANS  FROGRAN? 


THE  LAST  10  NVBC  VALUES i 


0.222000B 

02 

0.218000B 

02 

0.213000B 

02 

0.208000B 

02 

0.202000B 

02 

0.197000B 

02 

0.193000B 

02 

0.191000B 

02 

0.190000B 

02 

0.1B8000B 

02 

THE  25  PORECAST  VALUES: 

■ 

0.20389SB 

02 

0.200517B 

02 

0.197722B 

02 

0.195399B 

02 

0.19345SE 

02 

0.191S21B 

02 

0.190434B 

02 

0.1S924IB 

02 

0.188223B 

02 

0.107329B 

02 

0.18S541B 

02 

0.1S5B37B 

02 

0.185203B 

02 

0.184623B 

02 

0.1S4089B 

02 

0.183S91B 

02 

0.183121B 

02 

0.1E2S7SB 

02 

0.1B2249B 

02 

0.181E38E 

02 

0.181439B 

02 

0.181050B 

02 

0.180SS9B 

02 

0.180295B 

02 

0.179926B 

02 

THE  25  UPPER 

PORECAST  CONPIDBNCE  LIMITS: 

0.208094B 

02 

0.205056B 

02 

0.204779B 

02 

0.205033B 

02 

0.205462E 

02 

0.20S551B 

02 

0.207621B 

02 

0.208812B 

02 

0.210079B 

02 

0.211390E 

02 

0.212721E 

02 

0.214056B 

02 

0.21S381B 

02 

0.216687B 

02 

0.2179S9E 

02 

0.219221B 

02 

0.220442B 

02 

0.221830B 

02 

0.222783E 

02 

0.223903B 

02 

0.224989B 

02 

0.226043B 

02 

0.227064B 

02 

0.228055B 

02 

0;22901SB 

02 

ALPHA  POR  THE  CONPIDBNCE  LIMITS  XSt  0.100 


00  TOO  HAMT  TO  GO  AGAIN? 

_y 

ENTER  LETTER  POR  OPTION  YOU  WANT, 
c 

EXECUTION  BEGINS... 

SOLVES  POR  ROOTS  OP  CHARACTERISTIC  EQUATION. 

ENTER  NO.  OP  AR  PARAMETERS  (UNDIPP.  PORN). 

2 ' 

ENTER  AUTOREGRESSIVE  PARAMETERS: 

ENTER  PHI(l). 

1.8073 

ENTER  PHI (2). 

-0.8073 

THE  2 ROOTS  OP  CHARACTERISTIC  BQ.  AND  NORMS  ARE: 
1.239  0.0  1.239 

1.000  0.0  1.000 

00  YOU  NANT  TO  GO  AGAIN? 

y 


ENTER  LETTER  POR  OPTION  YOU  WANT, 
d 

EXECUTION  BEGINS... 

IS  YOUR  TINE  SERIES  SEASONAL? 
n 

ENTER  NUMBER  OP  NONSBASONAL  OIPPERBNCES. 

0 

00  YOU  MANT  TO  PLOT  AUTO  AND  PAUTO  OP  TRANSPORMBD  DATA? 
n 

OO  YOG  NANT  TO  GO  AGAIN? 

y 

ENTER  LETTER  POR  OPTION  YOU  WANT. 

s 

•EXECUTION  BEGINS... 

ENTER  NUMBER  OP  AR  PARAMETERS  (UNDIPPEREMCBD  PORM) 

2 

ENTER  ESTINATEO  AR  PARAMETER  PHI(l). 

1.8073 

ENTER  ESTIMATED  AR  PARAMETER  PHI (2). 

-0.8073 
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AR  PARAMETERS  ARB 
ARB  TRBSE  OK? 


ENTER  OVERALL  NA  CONSTANT 
>0. 00667 


enter  estimated  white  noise  VAR 
0.0177 

DO  YOU  WANT  starting  VALUES  TO  BE  TBS  LAST  VALUES 
OP  THE  ACTUAL  SERIES? 


START(l)-  19.000000 

START^(2)«  18.7999SS 

STARTIHC'VALUB(S)  OK? 

y 

ENTER  RANDOM  NUMBER  SEED  (BETWEEN  0 AND  1). 

0.12345 

ENTER  HUMBER  OP  VALUES  YOU  WANT  TO  SIMULATE  VIA  13 
SO 

DO  YOU  WANT  A PRINTOUT  OP  THE  PIRST  20  VALUES? 


SIMULATED  VALUES  ARE: 

18.50224  18.33957 
17.85420  17.69398 
16.94655  16.68689 
16.31497  16.64629 


17.98804 

17.26422 

16.26595 

17.62828 


18.16435 

17.50722 

li$;42499 

17.47137 


18.33694 

17.58267 

16.43227 

17.D4015 


DO  YOU  WANT  TO  GO  AGAIN? 


CONTROL  RBTORNEO  TO  CMS 
R»  - 


AUTOCORIiEUTlC^S  WITH  2 SIGMA  CANCS 
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AUTCCORREIAIICNS  UlTh  2 SIGM  BAKOS 
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APPENDIX  C 


KEYBOABD  SESSION  AND  OUTPUT  FROM  THE  ANALYSIS  OF 
TIME  SERIES  G (SEASONAL) 


This  appendix  contains  the  complete  terminal  keyboard 


session  emd  all  offline  output  from  the  analysis  of  Box 


and  Jenkins'  [Ref.  4]  time  series  G (International  airline 


passengers,  monthly  totals,  Jan  1949  - Dec  1960)  using  the 


Time  Series  Editor.  This  appendix  is  intended  to  supplement 


the  discussion  of  the  analysis  of  series  G given  in  Chapter 


VI  of  this  report. 


The  terminal  listings  are  presented  first,  followed  by 


the  computer  offline  output. 


repeat  l09in  nucepnt 

login  1621p44  400k 
ENTER  PASSMOROt 

ENTER  4-DIGIT  PROJECT  NUMBER  FOLLOWED  BY  4-CHARACTER  COST  CENTER  CODSt 
0444C172 

READY  AT  15.12.34  ON  09/16/78 
CMS  Version  ^.25 


link  2069p  191  192 
ENTER  PASSWORD! 

SET  TO  READ  ONLY 

ipl  CMS 

CMS  Version  3.25 

login  191 
R; 

login  192  t«p 
T (192)  R/0 
Rj 


tineser 

EACH  2 SECONDS  EXECUTION  TIME  IS  INDICATED  BY  *. 

YOU  HAVE  ENTERED  THE  TIME  SERIES  EDITOR.  r 

^ > • f p • ' ^ 

PLEASE  RESPOND  TO 'feXCH' OUERY 'WITH  AH  INPUT  AT  THE  TERMINAL. 

ENTER  ONLY  THE  FIRST  LSTTER  FOR  A WORD  RESPONSE. 

ENTER  NUMERICAL  VALUES  VIA  FORTRAN  FORMAT. 

TYPE  INTEGER  VALUES  (RIGHT  JUSTIFIED)"  FOR  NAMES  STARTING 

WITH  I THRU  H.  TYPE  FLOATING  VALUES.. WITH  DECIMAL  FOR  ALL  OTHERS. 

DO  YOU  WANT  A LIST  OF  THE  OPTIONS? 

y 

OPTION  DESCRIPTION 

GENERATE GENERATE  ANY  ARIMA  TIM^  ^BRIES 

AUTO CALCULATE  AUTOCORRELATIONS,  PAUTOS,  MEAN  AND  VARIANCE 

PLOT ‘-PLOT  A TIME  SERIES 

ESTIMATE  CALCULATE  MAX  LIKELIHOOD  ESTIMATES  OF  ARMA  PARAMETERS 

DIPF DIFFERENCE  A TI{IE  SERIES  • 

FORECAST  FORECAST  FUTUKE  VALUES,  CONSTRUCT  CONFIDENCE  INTERVALS 

TRANS  7-TRANSPORMS  VALUES  OF  A TIME  SERIES 

ROOTS  -DETERMINES  ROOTS  OF  ARIMA  CHARAfTTBRISTIC  EQUATION 

ZFORMAT ALTER  DATA  FILE  TO  FORMA*  6F15r< 

CMSWORK  PERFORM  CP/CMS  COMMANDS  IN  TIMESER  EXEC 

SIMULATE  SIMULATE  NONSEASONAL  TIME  SERIES 

YESTSEAS  -CALCULATE  INITIAL  SEASONAL  PARAMETERS 

WMARQRDT  MARQUARDT  SOLUTION  FOR  PARAMETER  ESTIMATES 

XSUMSQ  -CALCULATE  SUM  OF  SQUARES  FOR  ARBITRARY  PARAMETERS 

WOULD  YOU  LIKE  MORE  INFO?.  . - v'^' 

y ■ 

ENTER  OPTION  YOU  WANT  INFO  ABOUT. 

y 

YESTSEAS  THIS  PROGRAM  CALCULATES  INITIAL  PARAMETER  ESTIMATES 

FOR  SEASONAL  AND  NONSEASONAL  ARIMA  MODELS,  TO  BE  USED  AS  INPUTS  TO 
THE  WMARQRDT  PROGRAM. 


DO  YOU  WANT  INFO  ABOUT  ANOTHER  OPTION? 

y 

ENTER  OPTION  YOU  WANT  INFO  ABOUT. 
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MMARQROT  THIS  PROGRAM  CAI.CULATES  NON>bZNBAR  LEAST  SQUARES 

ESTIMATES  OP  BOX-JENKINS  PARAMETERS  FOR  A GENERAL  (SEASONAL 
OR  NONSBASONAL)  ARIMA  MODEL.  IT  REQUIRES  INITIAL  PARAMETER 
ESTIMATES  AS  STARTING  VALUES | THESE  MAY  BE  CALCULATED  USING 
P.tOGRAH  YBSTSEAS. 


DO  YOU  WANT  INFO  ABOUT  ANOTHER  OPTION? 

y 

ENTER  OPTION  YOU  WANT  INFO  ABOUT. 

X _ 

XSUNSQ  THIS  PROGRAM  CALCULATES  THE  SUM  OF  SQUARED  RESIDUALS 

FOR  ANY  SET  OF  ARINA  (SEASONAL  OR  NONSBASONAL)  PARAMETERS. 


DO  YOU  WANT  INFO  ABOUT  ANOTHER  OPTION? 
n 

*00  YOU  WANT  TO  TRY  A SESSION? 

y 

ENTER  LETTER  FOR  OPTION  YOU  WANT. 

C 

ENTER  DESIRED  CP/CHS  COMMANDS,  ONE  PER  LINE. 

WHEN  FINISHED  TYPBi  (GOTO  -QUES 
Stat 

P (I91)t  2«  FILES!  235  REC  IN  USB,  SI  LEFT  (oC  296),  79%  FULL  (2  CYL) 
cp  q £ 

FZLBSt>  NO  ROR,  NO  PRT,  NO  PUN 
listf  • et02eooi 

FILENAME  FILBTYPE  MODE  NO. REC.  DATE 
SERG  FT02F001  Pi  3 9/lS 

SERC  FT02F001  Pi  S • 9/lS 

alter  tecg  ftOlCOOl  pi  file  ftOlfOOl  pi 
(goto  -ques 

DO  YOU  WANT  TO  GO  AGAIN? 

y * 

ENTER  LETTER  FOR  OPTION  YOU  WANT. 

P 

IS  YOUR  DATA  IN  FILE  FT02F001? 

Y 

EXECUTION  BEGINS... 

ENTER  TITLE  FOR  PLOT 
secies  g data  / undlftecenced 
TIME  SERIES  PLOTS  HAVE  BEEN  PRINTED  OFFLINE 
DO  YOU  WANT  TO  GO  AGAIN? 

y 

ENTER  LETTER  FOR  OPTION  YOU  WANT, 
a 

IS  YOUR  DATA  IN  FILE  PT02F001? 

y 

EXECUTION  BEGINS... 

AUTOCORRELATIONS 


0.948 

0.87S 

0.807 

0.753 

0.714 

0.682 

0.663 

0.656 

0.671 

0.703 

0.743 

0.760 

0.713 

0.646 

0.586 

0.538 

0.500 

0.469 

0.450 

0.442 

0.457 

0.482 

0.517 

0.532 

0.494 

PARTIAL  AUTOCORRELATIONS 

0.948 

-0.229 

0.038 

0.094 

0.074 

0.008 

0.126 

0.090 

0.232 

0.166 

0.171 

-0.135 

-0.540 

-0.027 

0.091 

0.025 

0.032 

0.07: 

0.048 

-0.046 

0.04S 

-0.100 

0.052 

0.048 

-0.163 
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MEAN-  280.299 


VARIANCB  - 14292.0 


I ENTER  TITLE  FOR  PLOTS. 

autos  and  pautos  of  sac las  9 data  / undlffaranead 
TOOR  AOTO  ANO  PAOTO  PLOTS  RAVE  BEEN  PRINTED  OFFLINE. 
PICE  UP  IN  ROOM  X140  UNDER  YOUR  USER  10  NUHEBR. 
Cl  DO  YOU  WANT  TO  GO  AGAIN? 


ENTER  LETTER  FOR  OPTION  YOU  NANT. 

t 

IS  YOUR  DATA  IN  FILE  FT02F001? 

y ■ 

EXECUTION  BEGINS... 

DO  YOU  NANT^TO  TRANSLATE  THE  ORIGINi  N-X-SBIFT? 

n 

DO  YOU  MANT  TO  RESCALE  THE  VALUES t N-X*SCALB7 

n 

DO  YOU  WANT  A LOG  TRANSFORNATION7 

y 

transformation  is  H(I)-L0G(SCALE*(S(I)-SHIFT)*FACT0R)  WHBRBt 

SCALE*  1.00000  SHXFT-O.O  FACTOR«0.0 

YOUR  ORIGINAL  SERIES  IS  IN  DATA  FT02F001  PI. 

YOUR  TRANSFORMED  SERIES  IS  IN  FILE  FT02F001  PI. 


DO  YOU  WANT  TO  PLOT  THE  TRANSFORMED  VALUES? 

L ^ 

f EXECUTION  BEGINS... 

* ENTER  TITLE  FOR  PLOT 

plot  of  sacias  9 data  / natucal  I09  transfora 

TIME  SERIES  PLOTS  HAVE  BEEN  PRINTED  OFFLINE 

' DO  YOU  WANT  TO  GO  AGAIN? 

y 

I ENTER  LETTER  FOR  OPTION  YOU  WANT. 

i * 

' IS  YOUR  DATA  IN  FILE  FT02F001? 

y 

EXECUTION  BEGINS... 

anrocORRELATIONS 


1 0.954 

0.899  0.851  0.808 

0.779 

0.756 

0.738 

0.727 

0.734 

0.744 

1 1 0.758 

0.762  0.717  0.663 

0.618 

0.576 

0.544 

0.519 

0.501 

0.490 

i j 0.498 

0.506  0.517  0.520 

0.484 

i i 

1 1 

PARTIAL  AUTOCORRELATIONS 

' 0.954 

-0.118  0.054  0.024 

0.116 

0.044 

0.038 

0.100 

0.204 

0.064 

0.106 

-0.042  -0.485  -0.034 

0.042 

-0.044 

0.028 

0.037 

0.042 

0.014 

! 0.073 

-0.033  0.061  0.031 

-0.194 

• 

MEAN-  5.S4218  VARIANCE  ■ 0.193531 


! 


J 


ENTER  TITLE  FOR  PLOTS. 

i autos  and  pautos  of  sarias  9 data  / natural  I09  transfoca 

YOUR  AOTO  AND  PAUTO  PLOTS  HAVE  BEEN  PRINTED  OFFLINE. 

PICK  UP  IN  ROOM  1140  UNDER  YOOft  USER  ID  NUMBER. 

I DO  YOU  WANT  TO  GO  AGAIN? 

, y 
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t 





EMTmt  LETTBR  FOR  OPTION  YOU  HAHT. 
d 

IS  YOUR  DATA  IN  FILS  PT02F0017 

y 

EXECUTION  BEGINS... 

IS  YOUR  TIMS  SERIES  SEASONAi./ 

y 

ENTER  ORDER  OF  SEASONAL  DIFFERENCING. 

1 

ENTER  LENGTH  OF  SEASONAL  PERIOD  VIA  12. 

12 

ENTER  NUMBER  OF  NONSEASONAL  DIFFERENCES. 

1 

DO  YOU  WANT  TO  PLOT  AUTO  AND  PAUTO  OF  TRANSFORMED  DATA? 

y 

EXECUTION  BEGINS... 

AUTOCORRELATIONS 

-0.341  0.105  -0.202  0.021  0.05S  0.031  -O.OSS  -0.001  0.176  -0.076 

0.064  -0.387  0.152  -0.058  0.150  -0.139  0.071  0.016  -0.011  -0.117 

0.039  -0.091  0.223  -0.018  -0.100 

PARTIAL  AUTOCORRELATIONS 

-0.341  -0.013  -0.193  -0.125  0.033  0.035  -0.060  -0.020  0.226  0.043 

0.047  -0.339  -0.109  -0.077  -0.022  -0.140  0.026  0.115  -0.013  -0,167 

0.132  -0.072  0.143  -0.067  -0.103  , 


MBAN-0.290920B-03  VARIANCE  • 0. 2086048-02 


t i 


ENTER  TITLE  FOR  PLOTS. 

autos  and  pautos  for  sacias  g sOdata  / 1 ns  diff,  1 saas  dlff.  In  sfora 
YOUR  AUTO  AND  PAUTO  PLOTS  HAVE  BEEN  PRINTED  OFFLINE. 

PICK  OP  IN  ROOM  1140  UNDER  YOUR  USER  ID  NUMBER. 

DO  YOU  WANT  TO  GO  AGAIN7 

y 

ENTER  LETTER  FOR  OPTION  YOU  WANT, 
c 

ENTER  DESIRED  CP/CMS  COMMANDS,  ONE  PER  LINE. 

WHEN  FINISHED  TYPES  tCOTO  -CUES 
altar  fila  ftOlfOOl  pi  sava  ft03f001  pi 
altar  data  ft02f001  pi  sarg  ft02f001  pi 
stat 

P (191) t 29  FILES}  267  REC  IN  USB,  29  LEFT  (of  296),  90%  FULL  (2  CYL) 
aiasa  fila  ftOBfOOl 
cp  q £ 

FIl{ESt-  NO  RDR,  MO  PRT,  NO  PON 
(goto  -quas 

oa  YOU  WANT  TO  GO  AGAIN? 

y 

ENTER  LETTER  FOR  OPTION  YOU  HAHT. 

y 

IS  YOUR  DATA  IN  FILS  FT02F0017 

y 

EXECUTION  BEGINS... 

ENTER  LENGTH  OF  SEASON  VIA  FORMAT  12. 

12 

ENTER  NUMBER  OF  NON-SEASONAL  DIFFERENCES. 

1 

ENTER  NUMBER  OF  SEASONAL  DIFFERENCES. 

1 

ENTER  NUMBER  OF  NON-SEASONAL  AR  PARAMETERS. 

0 
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ENTER  NUMBER  OP  SEASONAL  AR  PARAMETERS 


ENTER  NUMBER  OP  MON-SEASONAL  MA  PARAMETERS 


ENTER  NUMBER  OP  SEASONAL  MA  PARAMETERS 


INITIAL  PARAMETER  ESTIMATES  POR  MARQRT  t 


0.390425 

0.534397 


THETA(l) 
THETAS (1) 


DO  YOU  WANT  TO  GO  AGAIN? 


ENTER  LETTER  POR  OPTION  YOU  MAMT 


IS  YOUR  DATA  IN  PILE  PT02P001? 


EXECUTION  BEGINS... 

ENTER  NUMBER  OP  NON-SEASONAL  DIPPERENCBS 


ENTER  NUMBER  OP  SEASONAL  DIPPERENCBS 


ENTER  NUMBER  OP  NON-SEASONAL  AR  PARAMETERS 


ENTER  NUMBER  OP  SEASONAL  AR  PARAMETERS 


ENTER  NUMBER  OP  NON-SEASONAL  MA  PARAMETERS 


ENTER  NUMBER  OP  SEASONAL  MA  PARAMETERS 


ENTER  LENGTH  OP  SEASON  VIA  FORMAT  12 
12 


NOW  INPUT  YOUR  INITIAL  PARAMETER  ESTIMATES,  AS  REQUESTED 


ENTER  NON-SEASONAL  MA  PARAMETER  THBTA(l) 
390425 

ENTER  SEASONAL  MA  PARAMETER  THETASd). 
534397 


CONVERGENCE  HAS  BEEN  REACHED  IN  MAX ( 7,  7)  ITERATIONS 


SELECTED  OUTPUT  POLLOWSt 


STANDARD  ERROR 


THBTA(l) 

THETASd) 


0.377152 

0.572387 


0.819933D-01 

0.780252D-01 


MOVING  AVERAGE  CONSTANT!  THETAO 


0.000291 


CHI-SQUARE  STATISTIC  POR  RESIDUAL  LACK  OP  PIT  - 29.713571 


DEGREES  OP  PREBOOM  - 38 


PROBABILITY  OP  EXCEEDING  STATISTIC  « 0.829438 


DO  YOU  NAMT  TO  PU)T  AOTO  AMO  PAOTO  OP  RESIDUALS? 

y 

EXECUTION  BEGINS... 

AUTOCORRELATIONS 


0.009 

0.026 

-0.129  -0.105 

0.078  0.077  -0.036  -0.033  0.103  -0.050 

0.026 

-0.020 

0.013  0.035 

0.067  -0.130  0.052  0.015  -0.093  -0.090 

-0.027 

-0.014 

0.213  0.009 

-0.042 

PARTIAL  AUTOCORRELATIONS 

0.009 

r 

0.026 

-0.130  -0.104 

0.088  0.068  -0.072  -0.029  0.151  -0.057 

-0.022 

0.018 

0.046  0.001 

0.053  -0.114  *.069  0.029  -0.129  -0.134 

0.042 

-0.014 

0.143  -0.025 

0.013 

HEAN-0.23600SE-02  VARIANCE  - 0.139713B*02 
ENTER  TITLE  POR  PLOTS. 

autos  and  pautos  oC  rasiduals  / (Of l,l)x(0, 1(1)12  nodal  of  sarias  9 
f'  YOUR  AOTO  AND  PAOTO  PLOTS  HAVE  BEEN  PRINTED  OPPLINE. 

PICK  UP  IN  ROOM  1140  UNDER  YOUR  USER  ID  NUMBER. 

* 

DO  YOU  WANT  TO  GO  AGAIN? 

^ y ^ . 

ENTER  LETTER  POR  OPTION  YOU  MANT. 

y 

IS  YOUR  DATA  IN  PILE  PT02P001? 

I y 

EXECUTION  BEGINS... 

I ENTER  LENGTH  OP  SEASON  VIA  PORMAT  12. 

12 

ENTER  NUMBER  OP  NON-SEASONAL  DIPPEREMCBS. 

1 

ENTER  NUMBER  OP  SEASONAL  DIPPERBNCES. 

1 

ENTER  NUMBER  OP  NON-SEASONAL  AR  PARAMETERS. 

1 

ENTER  NUMBER  OP  SEASONAL  AR  PARAMETERS. 

0 

ENTER  NUMBER  OP  NON-SEASONAL  NA  PARAMETERS. 

1 

ENTER  NUMBER  OP  SEASONAL  MA  PARAMETERS. 

1 


INITIAL  PARAMETER  ESTIMATES  POR  MARQRT  1 


I 


PHI(l)  « .0.112S99 
THETA(l)  ■ 0.490SS3 

THBTAS(l)  > 0.S33S79 

DO  YOU  WANT  TO  GO  AGAIN? 

y _ 

ENTER  LETTER  POR  OPTION  YOU  WANT, 
w 

IS  YOUR  DATA  IN  PILE  PT02P001? 

y 

EXECUTION  BEGINS... 

ENTER  NUMBER  OP  NON-SEASONAL  DIPPERENCES. 

1 

ENTER  NUMBER  OP  SEASONAL  DIPPERBNCES. 

1 

ENTER  NUMBER  OP  NON-SEASONAL  AR  PARAMETERS. 

1 


155 


ENTBR  NUMBER  OF  SEASONAL  AR  PARAMETERS. 

0 

ENTER  NUMBER  OP  NON-SBASONAL  MA  PARAMETERS. 

1 

ENTER  NUMBER  OP  SEASONAL  MA  PARAMETERS. 

1 

ENTER  LBNCTH  OP  SEASON  VIA  FORMAT  12. 

12 

NOW  INPUT  YOUR  INITIAL  PARAMETER  ESTIMATES,  AS  REQUESTED. 

ENTER  HON-SBASONAL  AR  PARAMETER  PHI(i). 

0.112699 

ENTER  NON-SBASONAL  MA  PARAMETER  THBTA(l). 

0.490883 

ENTER  SEASONAL  MA  PARAMETER  THBTAS(l). 

0.533679 

*** 

CONVERGENCE  BAS  BEEN  REACHED  IN  NAX(10,10)  ITERATIONS. 


SELECTED  OUTPUT  FOLLOWS t 


PARAMETER  ESTIMATE 


PHI(l)  - 0.146667 

THBTA(l)  ■ 0.509816 

TBBTAS(l)  > 0.573388 


STANDARD  ERROR 


0.2244640  00 
0.195934D  00 
0.7848960-01 


MOVING  AVERAGE  CONSTANT:  THETAI)  ■ 0.000248 


CHI-SQUARE  STATISTIC  FOR  RESIDUAL  LACK  OF  FIT  - 28.302543 
DEGREES  OF  FREEDOM  - 37 

PROBABILITY  OF  EXCEEDING  STATISTIC  ■ 0.847102 

DO  YOU  WANT  TO  PLOT  AUTO  AND  PAUTO  OF  RESIDUALS? 

y 

•EXECUTION  BEGINS... 

AUTOCORRELATIONS 

-0.013  0.073  -0.100  -0.088  0.077  0.073  -0.031  -0.033  O.IOS  -0.054 

0.033  -0.018  0.014  0.027  0.070  -0.133  0.051  0.008  -0.092  -0.084 

-0.020  -0.023  0.208  0.005  -0.036 

* 

PARTIAL  AUTOCORRELATIONS 

-0.013  0.073  -0.098  -0.096  0.092  0.081  -0.064  -0.040  0.152  -0.050 

-0.027  0.017  0.047  -0.004  0.053  -0.118  0.055  0.031  -0.123  -0.135 

0.051  -0.006  0.148  -0.019  0.003 


NEAN«0. 2530521-02 


VARIANCE  - 0.139373R-02 


ENTBR  TITLE  FOR  PLOTS. 

autos  and  pautos  of  casiduals  / (l,l,l)x(0, 1,1)12  aodal 
YOUR  AUTO  AND  PAUTO  PLOTS  HAVB  BEEN  PRINTED  OFFLINE. 
PICK  OF  IN  ROOM  1140  UNDER  YOOR  08RR  ID  NUMBER. 

DO  YOU  WANT  TO  GO  AGAIN? 

y 
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BMTBR  LBTTBR  fOR  OPTXOM  XOO  WANT. 

* . 

IS  YOUR  DATA  ZN  FZLB  rTOZrOOl? 

y 

EXECUTION  BEGINS.. . 

IS  YOUR  SERIES  SEASONAL? 

y 

ENTER  SEASONAL  DATA  AS  REQUESTED i 
ENTER  NUMBER  OF  SEASONAL  DIFFERENCES. 

1 

ENTER  NUMBER  OF  SEASONAL  AR  PARAMETERS. 

0 

ENTER  NUMBER  OP  SEASONAL  HA  FARAHBTBRS. 

1 

ENTER  LENGTH  OF  SEASON  VIA  FORMAT  Z2. 

12 

ENTER  NONSEASONAL  DATA  AS  REQUESTED i 
ENTER  NUMBER  OF  NON-SBASONAL  DIFFERENCES. 

1 

ENTER  NUMBER  OF  NON-SSASONAL  AR  PARAMETERS. 

0 

ENTER  NUMBER  OF  NON-SBASONAL  MA  PARAMETERS. 

1 

NOW  INPUT  foUR  PARAMETER  ESTIMATES > 


ENTER  MA  CONSTANT  TERM,  THETAO. 

.000291 

ENTER  NON-SEASONAL  MA  PARAMETER  THETA(l). 

.377152 

ENTER  SEASONAL  MA  PARAMETER  THETAS(l). 

.572387 

ENTER  MAXIMUM  FORECAST  LEAD  TIME  VIA  FORMAT  12. 

30 

ENTER  INDEX  FOR  PLOT  ORIGIN  VIA  FORMAT  13. 

100 

ENTER  INDEX  FOR  FORECAST  ORIGIN  VIA  FORMAT  13. 

131 

ENTER  SIGNIFICANCE  LEVEL  FOR  CONFIDENCE  INTERVALS. 

0.10 

WAS  YOUR  DATA  TRANSFORMED  ZN  THE  TRANS  PROGRAM? 

y 


DO  YOO  WANT  BASIC  OUTPUT  AT  THE  TERMINAL? 

y 


THE  LAST  10  WVEC  VALUES: 

0.418999E  03  0.4610018  03 

0.6060Q0B  03  0.507999E  03 

THE  30  FORECAST  VALUES: 


0.398S39E  03 
0.4697108  03 
0.460421B  03 
0.520593B  03 
0.70S666B  03 
0.526812E  03 


0.4I51B4E  03 
0.543456E  03 
0.405213E  03 
0.S07S48B  03 
0.590927B  03 
0.502002E  03 


0.472000B  03 
0.461001B  03 


0.395402B  03 
0.618194B  03 
0.448016R  03 
0.528791E  03 
0.5190898  03 
0.S87782B  03 


0.S3S001E  03 
0.390001B  03 


0.462697B  03 
0.6262768  03 
0.466864B  03 
0.611991B  03 
0.456978B  03 
0.S73220B  03 


0.622000B  03 
0.433002t  03 


0.45097ie  03 
0.524292E  03 
0.44474BB  'S 
0.696357E  03 
0.S05397B  03 
0.597386E  03 


157 


TUB  30  UPPER  FORECAST  CONFIDENCE  LIMITS: 


0.399S98E 

03 

0.416254E 

03 

0.396481E 

03 

0.463785E 

03 

0.452068B 

03 

0.470813B 

03 

0.544567B 

03 

0.619311E 

03 

0.627399B 

03 

0.525422E 

03 

0.461S57B 

03 

0.4063S4E 

03 

0.449173E 

03 

0.468030E 

03 

0.44S925B 

03 

0.S21779E 

03 

0.S08743E 

03 

0.529995E 

03 

0.613203B 

03 

0.697S77E 

03 

0.706894B 

03 

0.S92163E 

03 

0.520332E 

03 

0.458229B 

03 

0.506664B 

03 

0.528090B 

03 

0.S03292E 

03 

0.589083E 

03 

0.574S32B 

03 

0.S98708B 

0’ 

ALPHA  FOR  THE  COMFIOBNCB  LIMITS  IS:  0.100 
DO  TOO  HAMT  TO  GO  AGAIN? 

y 

ENTER  LETTER  FOR  OPTION  TOO  HAMT. 
f 

IS  TOUR  DATA  IN  FILE  FT02F0017 

y 

EXECUTION  BEGINS... 

IS  TOUR  SERIES  SEASONAL? 

y 

ENTER  SEASONAL  DATA  AS  REQUBSTEDt 

ENTER  NUMBER  OF  SEASONAL  DIFFERENCES. 

1 ' 

ENTER  NUMBER  OF  SEASONAL  AR  PARAMETERS. 

0 

ENTER  NUMBER  OF  SEASONAL  MA  PARAMETERS. 

1 

ENTER  LENGTH  OF  SEASON  VIA  FORMAT  12. 

12 

ENTER  NONSBASONAL  DATA  AS  REQUESTED: 

ENTER  NUMBER  OF  NON-SEASONAL  DIFFERENCES. 

1 

ENTER  NUMBER  OF  NON-SEASONAL  AR  PARAMETERS. 

1 

ENTER  NUMBER  OF  NON-SEASONAL  NA  PARAMETERS. 

1 

NOW  INPUT  YOUR  PARAMETER  ESTIMATES: 


ENTER  HA  CONSTANT  TERM,  THETAO. 

.000248 

ENTER  NON-SEASONAL  AR  PARAMETER  PNI(l). 

0.146667 

ENTER  NON-SEASONAL  MA  PARAMETER  THETA(l). 

0.S09B16 

ENTER  SEASONAL  MA  PARAMETER  THETAS (1). 

0.S73388 

ENTER  MAXIMUM  FORECAST  LEAD  TIME  VIA  FORMAT  12. 

30 

ENTER  INDEX  FOR  PLOT  ORIGIN  VIA  FORMAT  13. 

100 

ENTER  INDEX  FOR  FORECAST  ORIGIN  VIA  FORMAT  13. 

131 

ENTER  SIGNIFICANCE  LEVEL  FOR  CONFIDENCE  INTERVALS. 

0.10 

WAS  YOUR  DATA  TRANSFORMED  IN  THE  TRANS  PROGRAM? 

y 

DO  YOU  HAMT. BASIC  OUTPUT  AT  THE  TERMINAL? 

y 
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THE  LAST  10  WVEC  VALUES: 


f 


?■ 


j 

i 

{ 

[ 

j 

i 

f 

i 

i 

1 1 

.i 


1 ‘ 


0.418999E 

03 

0.461001B 

03 

0.472000E 

03 

0.535001B 

03 

0.822000E 

03 

0.606000B 

03 

0.S07999E 

03 

0.461001B 

03 

0.390001C 

03 

0.432002B 

03 

THE  30  FORECAST  VALUES: 

0.398650E.  03 

0.414875E 

03 

0.395125B 

03 

0.462366B 

03 

0.450651E 

03 

0.469348B 

03 

0.543045E 

03 

0.617700E 

03 

0.625753E 

03 

0.523897E 

03 

0.4S0071B 

03 

0.404899B 

03 

0.447739B 

03 

0.466360E 

03 

0.444325E 

03 

0.520098B 

03 

0.507066B 

03 

0.528257B 

03 

0.611381C 

03 

0.695632E 

03 

0.704905B 

03 

0.59033SB 

03 

0.518566E 

03 

0.456512E 

03 

0.504960B 

03 

0.526113B 

03 

0.501400E 

03 

0.587074B 

03 

0.572532E 

03 

0.596632B 

03 

THE  30  UPPER 

FORECAST  CONFIDENCE  LIMITS: 

0.399709B 

03 

0.415945B 

03 

0.398204B 

03 

0.463452B 

03 

0.451744E 

03 

0.470448B 

03 

0.544151B 

03 

0.618812B 

03 

0.626870B 

03 

0.525020B 

03 

0.4S1199B 

03 

0.406032E 

03 

0.448887E 

03 

0.4875178 

03 

0.445491C 

03 

0.S21270E 

03 

0.508249B 

03 

0.5294488 

03 

0.812579B 

03 

0.696837B 

03 

0.706117E 

03 

0.S91S56E 

03 

0.519792B 

03 

0.45774SB 

03 

0.5082078 

03 

0.527371B 

03 

0.502668B 

03 

0.588353E 

03 

0.573821B 

03 

0.597930B 

03 

ALPHA  POE  THE  CONPXOBNCB  LIMITS  IS:  0.100 
00  you  WANT  TO  GO  AGAIN? 

y 


ENTBR  LBTTBR  FOR  OPTION  YOU  WANT. 
f 

IS  YOUR  DATA  IN  PILE  PT02F0017 

y 

EXECUTION  BEGINS... 

IS  YOUR  SERIES  SEASONAL? 

y 

ENTER  SEASONAL  DATA  AS  RBQUBSTBD: 

ENTER  NUNBBR  OP  SEASONAL  OIPPBRBNCES. 

1 

ENTER  NUMBER  OP  SEASONAL  AR  PARAMETERS. 

0 

ENTER  NUMBER  OP  SEASONAL  NA  PARAMETERS. 

1 

ENTER  LENGTH  OP  SEASON  VIA  FORMAT  12. 

12 

ENTER  NONSEASONAL  DATA  AS  REQUESTED t 
ENTER  NUMBER  OP  NON-SEASONAL  0IPPERBNCB8. 

1 

ENTER  NUMBER  OP  NON-SEASONAL  AR  PARAMETERS. 

0 

ENTER  NUMBER  OF  NON-SEASONAL  NA  PARAMETERS. 

1 


NON  INPUT  YOUR  PARAMETER  ESTIMATES: 


ENTER  MA  CONSTANT  TERM,  THETAO. 

.000291 

CUTER  NON-SEASONAL  HA  PARAMETER  THETA(l). 
.377152 

CUTER  SEASONAL  NA  PARAMETER  THETAS(l). 

.572387 

ENTER  MAXIMUM  FORECAST  LEAD  TINE  VIA  FORMAT  12. 

40 
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ENTER  INDEX  FOR  PLOT  ORIGIN  VIA  FORMAT  13. 

075 

ENTER  INDEX  FOR  FORECAST  ORIGIN  VIA  FORMAT  13. 

120 

ENTER  SIGNIFICANCE  LEVEL  FOR  CONFIDENCE  INTERVALS 

0.10 


WAS  TOUR  DATA  TRANSFORMED  IN  TBB  TRANS  PROGRAM? 


DO  YOU  MANT  BASIC  OUTPUT  AT  THE  TERNINAL? 


ENTER  LETTER  FOR  OPTION  YOU  WANT 


IS  YOUR  DATA  IN  FILE  FT02F001? 


EXECUTION  BEGINS... 

IS  YOUR  TIME  SERIES  SEASONAL? 


ENTER  ORDER  OF  SEASONAL  DIFFERENCING 


ENTER  LENGTH  OF  SEASONAL  PERIOD  VIA  12 


ENTER  NUMBER  OF  NONSEASONAL  DIFFERENCES 


DO  YOU  WANT  TO  PLOT  AUTO  AND  PAUTO  OF  TRANSFORMED  DATA? 


DO  YOU  WANT  TO  GO  AGAIN? 


ENTER  LETTER  FOR  OPTION  YOU  WANT 


ENTER  DESIRED  CP/CMS  COMMANDS,  ONE  PER  LINE 
WHEN  FINISHED  TYPEt  (GOTO  -CUES 
alter  file  £t02C001  pi  Inserg  £t02f001  pi 
alter  file  £t03f001  pi  file  £t02£001  pi 
tgoto  -ques 

DO  YOU  WANT  TO  GO  AGAIN? 


ENTER  LETTER  FOR  OPTION  YOU  WANT 


IS  YOUR  DATA  IN  FILE  FT02F001? 


EXECUTION  BEGINS 


IS  YOUR  SERIES  SEASONAL? 


ENTER  LENGTH  OF  SEASON  VIA  FORMAT  12 


ENTER  NUMBER  OF  SEASONAL  AR  PARAMETERS 


ENTER  NUMBER  OF  SEASONAL  MA  PARAMETERS 


ENTER  NUMBER  OF  NON-SEASONAL  AR  PARAMETERS 


ENTER  NUMBER  OF  NON-SEASONAL  MA  PARAMETERS 


NOW  INPUT  YOUR  INITIAL  PARAMETER  ESTIMATES.  AS  REQUESTED 


ENTER  NON-SEASONAL  MA  PARAMETER  THBTA(l) 
377152 

ENTER  SEASONAL  MA  PARAMETER  THETAS (1). 
5723S7 


FOR  THB  GIVEN  INPUT  PARAMETERS  AND  MODEL 

THE  RESIDUAL  SOM  OP  SQUARES  ISi  0.1S192935D  00 


DO  YOU  WANT  TO  TEST  DIFFERENT  PARAMETER  VALUES? 


DO  YOU  WANT  TO  TEST  A DIFFERENT  MODEL? 


IS  YOUR  SERIES  SEASONAL? 

ENTER  LENGTH  OF  SEASON  VIA  FORMAT  12. 

12 

ENTER  NUMBER  OF  SEASONAL  AR  PARAMETERS. 

) 

ENTER  NUMBER  OF  SEASONAL  MA  PARAMETERS. 

I 

ENTER  NUMBER  OF  NON-SEASONAL  AR  PARAMETERS. 

I 

ENTER  NUMBER  OF  NON-SEASONAL  MA  PARAMETERS. 

I 

NOW  INPUT  YOUR  INITIAL  PARAMETER  ESTIMATES , AS  REQUESTED. 

ENTER  NON'SEASONAL  AR  PARAMETER  PHI(l). 

.146667 

ENTER  NON'SEASONAL  MA  PARAMETER- TaBTA(l) . 

.509S16 

ENTER  SEASONAL  MA  PARAMETER  THETAS <1). 

.573388 


FOR  THE  GIVEN  INPUT  PARAMETERS  AND  MODEL 

THB  RESIDUAL  SUM  OF  SQUARES  ISi  0.18164378D  00 


DO  YOU  WANT  TO  TEST  DIFFERENT  PARAMETER  VALUES? 
n 

DO  YOU  WANT  TO  TEST  A DIFFERENT  MODEL? 
n 

DO  YOU  WANT  TO  GO  AGAINt 
n 

CONTROL  RETURMBO  TO  CMS 


X'SCAlEt  0.179E  01  UNITS 

Y-SC^LEt  "**•  0.863E  01  UNITS 
SeHeS  6 C*T*  / UNOIFPERENCEO 
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AVTCCORREIATIONS  M|T|-  2 SlCf*  BANOS 


mOCORRI  LA  TICKS  WITH  2 SICf>A  BANDS 


Ua'.'i'-V'H  tL'IVK  l'JlbULl->U;VJU:UlU^UMl.'lbU:W  lU 

iMr«<»ii»>w«nn(weD>4tf>f»Orw''ir>p><or>o>^«iniwitt 

eg^«*<o«>‘»o*»op‘oog»oo>«»»CTC>w*<^v<«'«w*><«* 

r»eo<^o^rtOooooso-jooe‘Oo<^^^ 


III  I I I I I I I I I I I I I I I I I I I I I 

••••••••••••••••««•••••••• 

O^JO^OOOOOOOf^OOoOoOOOOOOOOO 
II  II  I III  III  I 


Mc^o^oooooo 
I I T I • I I • I I I 


I I I I I I I I I I I I I I I t I • I I • I • I 

ui;s»*u'u*«ijui«»  tt:;;iti;u;iV|i-U4«^uai.‘u^aixujik.uAutj 

ooo^99^o«9»o^0r;9»09(>> 

o of  mtr>  ^ f**  O' o ^ ^ to  oo a>  c>  oi  Of  <^wo  ^ <M 

1*40*1000 

iM^-tOfMO»*r<tv-*oi*r<fi>^*4^^f#^'rrn<no<»o-.Ooi^**4 


p >Oii  *P*<|  t,  I 

0 or^ooooooo^oooowoo-oooowoo 

1 I I I I I t I I I I I I I I I I I • I I I I 
u:  U*Wtt;i4fUL*M*l«Mbl*/IUi«  ijatUiUitdiOM 

4%  r»ou%OP^r*«f«ii4r^«D^rCO'<fco4MO(^«nr**0*i4<o 
m <oio*iPOaKJ^'0«^#«^^0'0;ov*9*40*.^V*^o 

tf*  CM*^«-«9(^(nr400fW(M09iOtf«C«^^On*4l9<M<0 
10  o*roo>’-M.-4>0'Ow*««o»mO'no'*>*^'r<n4>ooio 
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I . APPENDIX  D 

PROGRAM  LISTINGS 

This  Appendix  contains  progr2un  listings  of  all  progreuns 
resident  in  the  Time  Series  Editor.  It  does  not  include 
listings  of  subroutines  used  from  other  libraries,  such 
as  IMSL  routines. 
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APPENDIX  E 


LISTING  OF  DATA  SETS  RESIDING 
IN  THE  TINE  SERIES  EDITOR 


This  appendix  contains  listings  of  the  data  sets  now 
residing  in  the  Time  Series  Editor.  The  data  consists  of 


Box  and  Jenkins'  [Ref.  4]  time  series  C (Ref.  4,  p.  528] 


and  time  series  G [Ref.  4,  p.  531],  as  well  as  a data  set 


containing  monthly  Monterey ^ California,  rainfall  data 


from  January  1931  through  December  1976  . 
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